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Abstract 



The perturbative calculation of the lifetime of fermion excitations in a QED 
plasma at high temperature is plagued with infrared divergences which are not 
eliminated by the screening corrections. The physical processes responsible for these 
divergences are the collisions involving the exchange of longwavelength, quasistatic, 
magnetic photons, which are not screened by plasma effects. The leading divergences 
can be resummed in a non-perturbative treatement based on a generalization of 
the Bloch-Nordsieck model at finite temperature. The resulting expression of the 
fermion propagator is free of infrared problems, and exhibits a non- exponential 
damping at large times: Sji{t) ~ exp{— aTt Inuipt}, where cOp = eT/3 is the plasma 
frequency and a = e^/47r. 
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1 Introduction 



The study of the elementary excitations of ultrarelativistic plasmas, such as the quark- 
gluon plasma, has received much attention in the recent past [1 — 11] (See also ||12|, |l^ 
for recent reviews and more references.) The physical picture which emerges is that of a 
system with two types of degrees of freedom: i) the plasma quasiparticles, whose energy 
is of the order of the temperature T; ii) the collective excitations, whose typical energy 
is gT, where g is the gauge coupling, assumed to be small: (7 <^ 1 (in QED, g = e is the 
electric charge). For this picture to make sense, however, it is important that the lifetime 
of the excitations be large compared to the typical period of the modes. 

Information about the lifetime is obtained from the retarded propagator. A usual 
expectation is that Sji{t, p) decays exponentially in time, Sji{t, p) ~ Q-^E{p)t^-'y{p)t ^ th.dX 
l*S'ij(t,p)p ~ e~'"*^^^* with T{p) = 2'j{p), which identifies the lifetime of the single particle 
excitation as r(p) = l/r{p). The exponential decay may then be associated to a pole of 
the Fourier transform Sr{lj,p), located at u; = E{p) — i'y{p). The quasiparticles are well 
defined if their lifetime r is much larger than the period ~ 1/i? of the field oscillations, 
that is, if the damping rate 7 is small compared to the energy E. If this is the case, 
the respective damping rates can be computed from the imaginary part of the on-shell 
self-energy, T,{uj = E{p),p). Such calculations suggest that 7 ~ g'^T [|, ^ for both the 
single-particle and the collective excitations. In the weak coupling regime g <^ 1, this is 
indeed small compared to the corresponding energies (of order T and gT, respectively), 
suggesting that the quasiparticles are well defined, and the collective modes are weakly 
damped. However, the computation of 7 in perturbation theory is plagued with infrared 
divergences, which casts doubt on the validity of these statements 0, [15 — 26] 

The first attempts to calculate the damping rates were made in the early 80 's. It 
was then found that, to one-loop order, the damping rate of the soft collective excita- 
tions in the hot QCD plasma was gauge-dependent, and could turn out negative in some 
gauges (see Ref. [^] for a survey of this problem). Decisive progress on this problem was 
made by Pisarski and Braaten and Pisarski who identified the resummation needed 
to obtain the screening corrections in a gauge-invariant way 0] (the resummation of the 
so called "hard thermal loops" (HTL)). Such screening corrections are sufficient to make 
finite the transport cross-sections 0, and also the damping rates of excitations with 
zero momentum]^, ^. At the same time, however, it has been remarked^ that the HTL 
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resummation is not sufficient to render finite the damping rates of excitations witli non 
vanisfiing momenta. Tlie remaining infrared divergences are due to collisions involving 
the exchange of longwavelength, quasistatic, magnetic photons (or gluons), which are not 
screened in the hard thermal loop approximation. Such divergences affect the computa- 
tion of the damping rates of charged excitations, in both abelian and non-abelian gauge 
theories. Thus, in the lowest order calculations of Refs. [|], [15 — 26], one meets the same 
logarithmic divergence for electrons in QED, for charged scalars in SQED, and for quarks 
and gluons in QCD. (There is no such problem for the photon damping rate, which is IR 
finite and of order g'^T since photons do not couple directly to gluons or to them- 
selves.) Furthermore, the problem appears for both soft (p ~ gT) and hard (p ~ T) 
quasiparticles. In QCD this problem is generally avoided by the ad-hoc introduction of an 
IR cut-off ( "magnetic screening mass" ) ~ g'^T, which is expected to appear dynamically 
from gluon self-interactions E^. In QED, on the other hand, it is known that no magnetic 



screening can occur ^0[, so that the solution of the problem must lie somewhere else. 

In order to make the damping rate 7 finite, Lebedev and Smilga proposed a self- 
consistent computation of the damping rate 7|]T3[, by including 7 also in internal propaga- 
tors. However, the resulting self-energy is not analytic near the complex mass-shell, and 
the logarithmic divergence actually reappears when the discontinuity of the self-energy is 
evaluated aX uj = E — i'-y [|l^, |18[. More thorough resummations of the fermion line led 
to the conclusion that the full fermion propagator has actually no quasiparticle pole in 



the complex energy plane [p3|, These analyses left unanswered, however, the question 
of the large time behavior of the retarded propagator. As we have shown in a previous 
letter the answer to this question requires resummations for both the fermion prop- 
agator and the photon-electron vertex function. Such resummations modify the analytic 
structure of the retarded propagator: indeed, as we shall see, they make it analytic in the 
vicinity of the mass-shell. 

The need for a nonperturbative analysis follows from the fact that infrared diver- 
gences occur in all orders of perturbation theory. The leading divergences arise, in all 
orders, from the same kinematical regime as in the one loop calculation, namely from the 
exchange of soft quasistatic magnetic photons. In the imaginary time formalism, these 
divergences are concentrated in diagrams in which the photon lines carry zero Matsubara 
frequency (to be referred as static modes in what follows). In this sense, they appear as 
the divergences of an effective three-dimensional gauge theory, which is intrinsically non- 
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perturbative. Still, this effective "dimensional reduction" brings in simplifications which 
can be exploited to arrive at an explicit solution of the problem. 

We concentrate in this paper on the damping rate of fermionic excitations in hot 
QED plasmas. Our analysis is based on the Bloch-Nordsieck (or eikonal) approximation p2 
At zero temperature, this approximation provides an all-order solution to the infrared 
catastrophe, and correctly describes the mass-shell structure of the 4-dimensional fermion 
propagator ||33|| . At finite temperature, the Bloch-Nordsieck approximation has been pre- 
viously used, by Weldon, to verify the cancellation of the infrared divergences in the 
production rate for soft real photons |3^. Let us also mention that an attempt to solve 
the IR problem of the damping rate, using the BN approximation in the same spirit as 
in the present paper, has been reported in Ref. [BH]. However, although the final re- 



sult obtained in |35| is similar to ours, the derivations there are plagued with several 



inconsistencies, some of which are pointed out in |31 



In this paper, we shall consider (in section 3) a different generalization of the Bloch- 
Nordsieck (BN) model at finite temperature, which is better suited to study the infrared 
structure of the fermion propagator. Our approach is a natural extension of the method 
used in Ref. iSH] in QEDg^^ at zero temperature. However, the resulting imaginary-time 
BN propagator does not exponentiate in an obvious way, and thus cannot be written in 
closed form, in contrast to the usual, zero-temperature propagator. Still, we can obtain an 
explicit solution once we restrict ourselves to the static Matsubara photon modes. We thus 
get the retarded propagator Sii(t,p), and study its large time behavior (section 4). The 
final result is that, for times t ^ 1/ gT, the propagator does not show the usual exponential 
decay alluded to before, but the more complicated behavior Sn^t, p) ~ e~'-^(p)*e~""^* inujpt^ 
where Up ~ gT is the plasma frequency, and i?(p) ~ p ~ T is the average energy of 
the hard fermion. This corresponds to a typical lifetime ~ g'^Tln{l/g), which is 
similar to the one provided by the perturbation theory with an IR cut-off of the order 
g'^T. Since, as t oo, S'/j(t) is decreasing faster than any exponential, the Fourier 
transform of Sfi{t,p), Sii{u!,p), is an entire function in the complex energy plane. The 
existence of the quasiparticle is therefore not signaled by the presence of a pole of Sr^u) 
in the complex energy plane. However, the associated spectral density has the shape 
of a resonance strongly peaked around u = E{p), with a typical width of the order 
1/r ~ g'^T ln{l / g) . With minor modifications, the above conclusions also apply for the 
soft (collective) excitations, with momenta p ~ gT, whose lifetimes are found to depend 
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on the group velocities |f±| < 1 (section 5). 

At this stage it is useful to specify the notations and the conventions to be used 
throughout. The analytic propagator is defined in the complex energy plane by the 
spectral representation 

^(-,p) = r^pjp^. (1.1) 

The Matsubara propagator is obtained from eq. ( |1.1|) by setting u = iun, with Un = 
(2n+l)7rT and integer n. At tree level, pf(p°, p) = ^Po{p^,p) where = p'^'jfj,, €p = |p| = p, 
and 

Poip',p) = ^ {Sip' - e,) - 5{p' + e,)), (1.2) 

so that 



— — = P + — — 

UJ"^ — p^ uo — p 00 + p 



So{u, p) = — 2 _ 2 = ~ — I ^+(P) + h-ip), (1-3) 



where h±{'p) = (7° =F P ■ 7)/2, with p = p/p. 

The full fermion propagator is given by the Dyson-Schwinger equation 

S-\u;,p) = So\u;,p) + S(^,p). (1.4) 

The most general form of the self-energy S which is compatible with the rotational and 
chiral symmetries is 

S(tu,p) = a(a;,p)7° + 6(u;,p)p ■ 7 = /i_(p)S+(cj,p) - /i+(p)S_(u;,p), (1.5) 

where 

S±(a;,p) = ±^tr(/i±(p)E(a;,p)). (1.6) 
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Using this decomposition of S onto h±, and the analogous one for 5*0, eq. ( |1.3| ), one can 
easily invert eq. ( |1.4| ) to get the full propagator: 

S{uj,p) = A+{uj,p)h4p) + A4u,p)h4p), (1.7) 

where 

A±{uj,p) 



io =F {p+Y:±{uj,p)) 
4 




P k 

Figure 1: The one- loop fermion self-energy 

The retarded propagator is obtained as the boundary value of the analytic prop- 
agator ( |1.1|) when oj approaches the real axis from above, i.e., Sr{uj^^) = S{uj + iri,p), 
where u is real and ^ 0+. In the time representation, 

/°° rlij) r°° r\(j) 

^e-^-'Sn{u;,p) = i9{t) / Vf(^,P)- (1-9) 

-oo zvr J-oo zvr 

The large time behavior of S'ij(t, p) is determined by the analytic structure of Sr{uj,p) 
when continued to complex values of uj. In the upper half plane, Sr{ijj) coincides with 
the analytic propagator (|1.1D. In the lower half plane, Sr{(jj) is defined by continuation 
across the real axis, and it may have there singularities. The large time behavior of 
Suit) is controlled in most cases by the singularity of Sr{uj) which lies closest to the 
real axis. If this is located aX uj = E{p) — i^{p), then ^^^(t, p) ~ /(t, p) e~*^'^^^*e"'''^^^*, 
where the prefactor /(t, p) is slowly varying, and depends on the specific nature of the 
singularity. This conventional picture breaks down in gauge theories since, as we shall 
discuss in the next section, the perturbative estimate of 7 turns out to be IR divergent. 
The resummation of the leading infrared divergences, carried out in section 3, produces 
a propagator which has no singularity in the complex u plane. We shall then find it 



convenient to calculate Sji{t) directly, rather than from the Fourier transform (|1.9|). 



2 The one-loop damping rate for the hard fermion 

In this section, we review the perturbative calculations of the damping rate for a hard 
fermion, with momentum p ~ T [15 — 26] focussing on the infrared divergences which 
arise in such calculations. We assume here, as customary, that the dominant singularity 
of the retarded propagator is a simple pole whose location goes back into the tree-level 
pole at u = p when g 0. 
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p p-q 

Figure 2: Two- loop diagram contributing to the fermion self-energy 
2.1 Physical interpretation of the damping 

To leading order in g, the self-energy is given by the one-loop diagram in Fig. |I|. This gives 
no contribution to the damping rate 7. Indeed, when evaluated on the free mass-shell, i.e. 
at uj = p, the imaginary part of the one-loop self-energy vanishes because of kinematics. 
(At finite temperature this argument involves subtleties which are discussed in Appendix 
B.) 

The leading contribution to 7 comes therefore from the two-loop diagram in Fig. 
and turns out to be quadratically infrared divergent (see, e.g., Refs. [^, ^ |T^, [T^). 



The on-shell imaginary part is obtained by cutting the diagram in Fig. |^ through the 
internal fermion loop and the lower fermion propagator. Physically, this imaginary part 
accounts for the scattering of the incoming electron (with four momentum = (ep, p) and 
Cp = p) off a thermal fermion (electron or positron), calculated in the Born approximation 
(see Fig. 3). The total interaction rate is given by 

1 /■ 

r(p) = ^ y dpi dp2 dp3 {271^ 6'-'^\p + pi-p2~ P3) 

{n,{l - n2)(l - ria) + (1 - ^1)722^3} \M\^ (2.1) 

and coincides with twice the damping rate 7(p), as computed from the two-loop self- 
energy in Fig. 2: T{p) = 2'~f{p). This identity extends to finite temperature the usual 
physical interpretation of the self-energy discontinuity in terms of cross-sections for phys- 
ical processes, and can be verified through an explicit calculation (see also below). 
The notations in eq. ( p.l|) are as follows: all the particles are on the mass-shell (i.e., e = p 
and ej = pi for i = 1, 2, 3), and we have denoted / dpi = J{d^pi/{2n)^2ei). The factors 
Hi = n{ei) are the thermal occupation numbers for fermions (n(e) = l/(e^'^ -|- 1)). Note 
that, for fermions, the rates of the direct and of the reverse processes have to be added to 
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Figure 3: Fermion-fermion elastic scattering in the Born approximation 



give the total depopulation of the fermion state with momentum p'^ |T^. Finally, is 
the scattering matrix element squared, averaged over the spin s of the incoming electron, 
and summed over the spins si, S2 and S3 of the other three particles. In the Born approx- 
imation, |A^p is independent of the temperature and involves the propagator D^y{q) of 
the exchanged photon (with q = p — P3 = P2 — Pi) ■ Specifically P5|, 



\M\' = lQg'D^MD;M{p'p', + P^f - g^'iP ■ P3)}{PiP2 + P2P1 - 9'\Pi ■P2)}. (2.2) 

We shall use below the Coulomb gauge where the only non-trivial components of D^u{q) 
are the electric (or longitudinal) one Doo(q') = and the magnetic (or transverse) 

one Dij{q) = {6ij - qiqj)At{q). 

Since the interaction rate ( p.lD is dominated by soft momentum transfers q <^ 
while the external momenta are typically of the order of T, we can simplify the matrix 
element |A^p by setting p ~ ps and pi ^ p2 in eq. (|2.2|) , and obtain 



|Alp ~ Ug^pl |Az(g) + (v x q) ■ (vi x q)Ai(g)|', (2.3) 

with V = p and Vi = pi. Furthermore, we use energy conservation to write go = ^ ~ ^3 = 
62 — 61, that is, 

9o = p - |p - q| = |pi + q| - Pi , 

which, for g <^ T, becomes 

go~v-q~Vi-q. (2.4) 
The statistical factors in eq. ( ^.1[ ) satisfy the following identity: 

ni(l - n2){l - ns) + (1 - ni)n2n3 = {ui - n2)(l + A^(go) - %) , (2.5) 



which features iV(go), the Bose- Einstein thermal factor for the virtual photon. Since 
£2 = ei + go, and go < Pi ~ T, 

(\ti (\ti 
[n,-n,){l + N{q,)-n,) ^ -—q,N{q,) ^ -T— , (2.6) 

where we have used the fact that, at small go ^ 

l + iV(go)-n3 ~ iV(go) ~ T/q,. (2.7) 

Finally, we use eq. ( |2.4| ) to rewrite the integrations over p2 and pa as follows 

d^P2 M^P3 (2^)4^(4) (^ + ^^_^^_ 



(27r)3 7 (27r)3 

SO that we may use pi, q and go as independent integration variables in eq. (^.l]): 

5(go - V ■ q) 5(go - Vi ■ q) |A;(g) + (v x q) ■ (vi X q)Ai(g)|'. (2.9) 

We perform the angular integrations over vi = pi and q by using the delta functions, 
while the radial integration over pi gives 

r n(dn\ TT^T^ 



We obtain finally 



^1-3 |A*(go,g)r , (2.11) 



6 Jo J-q 27r 2 Y g 

where the upper cut-off g* distinguishes between soft and hard momenta: gT ^ g* ^ T. 
Since the g-integral is dominated by IR momenta, its leading order value is actually 
independent of g*. 

The two terms within the parentheses in eq. ( |2.11| ) correspond to the exchange of an 
electric and of a magnetic photon respectively. For a bare photon, we have |A/(go, g)p = 
1/g^ and |Ai(go, g)p = l/(go — g^)^, so that the g-integral in eq. (|2.11|) shows a quadratic 
IR divergence: 

An Jo g3 



r-^/;^. (2.12) 
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This divergence reflects the singular behaviour of the Rutherford cross-section for forward 



scattering p8[] . 

As well known, however, the quadratic divergence is removed by the screening cor- 
rections contained in the photon polarization tensor. These modify the electric and mag- 
netic propagators as follows 

*Ai{qo,q) = — r, *At{qo,q) = — ^ r, (2.13) 

q^-6Ui{qo,q) - - S^t{qo, q) 

where SUi and 6Ilt are the respective pieces of the photon polarisation tensor (in the hard 
thermal loop approximation [|l], |^). We shall see below that the leading IR contribution 
comes from the domain qQ <^ q <^T, where we can use the approximate expressions 

SUiiqo « q) ~ 3ul = ml, 6Ut{qo <^ q) - -i ^ ujI^ . (2.14) 
f A ' q 

We see that screening occurs in different ways in the electric and the magnetic sectors. In 
the electric sector, the familiar static Debye screening provides an IR cut-off rrio ^ gT. 
Accordingly, the electric contribution to F is finite, and of the order Ti ~ g'^T^/ml ~ g^T. 
Its exact value can be computed by numerical integration |T^. In the magnetic sector, 
screening occurs only for nonzero frequency go [@) 0] • This comes from the imaginary part 
of the polarisation tensor, and can be associated to the Landau damping of space-like 
photons (go < q^)- This "dynamical screening" is not sufficient to completely remove the 
IR divergence of Vt : 

g^T^ r* r dgo 1 



12 Jo J-q 2ti g4 + (37rcj2go/4g)2 

^'■•^'a.ctanf^U'!! r^!?. (2.15) 



ir'^ Jo q \ Aq^ J 2t[ Jo q 

In writing the last equality, we payed attention only to the dominant, logarithmically 
divergent, contribution. To isolate it, we have written 



arctan 



p 



" 2 ' 



as appropriate for g <S Up, and we have introduced the upper cut-off Up ~ gT to approx- 
imately account for the correct UV behaviour of the integrand: namely, as g ^ oOp, the 
integrand is decreasing like uj^/q^, so that the g-integral is indeed cut-off at g ~ Up. 

The remaining IR divergence in eq. ( |2.15|) is due to collisions involving the exchange 
of very soft (|q ^ 0|), quasistatic (go 0) magnetic photons, which are not screened by 
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plasma effects. To see that, note that the IR contribution to Tt comes from momenta 

q <^ gT, where \At{qQ, g)p is almost a delta function of Qq: 

1 4 
|At(go,g)|^ - — — ^g^o ^ — ^^(go)- (2.16) 

This is so because, as go 0, the imaginary part of the polarisation tensor vanishes 
linearly (see the second equation ( |2.14| )), a property which can be related to the behaviour 
of the phase space for the Landau damping processes. Since energy conservation requires 
qo = q cos 6, where 6 is the angle between the momentum of the virtual photon (q) and 
that of the incoming fermion (p), the magnetic photons which are responsible for the 
singularity are emitted, or absorbed, at nearly 90 degrees. 

To conclude this subsection, we note that, if we temporarily leave aside the loga- 
rithmic divergence, then both the electric and the magnetic damping rates are of order 
g'^T, rather than g^T as one would naively expect by looking at the diagrams in Figs. 2 
and 3. This situation has been sometimes referred as anomalous damping [|15 , and is a 
consequence of the strong sensitivity of the scattering cross section to the IR behavior 
of the photon propagator. By comparison, the other two-body collisions leading to the 
damping of the fermion, namely the Compton scattering and the annihilation process, 
are less IR singular — as they involve the exchange of a virtual fermion — and only 
contribute at order g'^T. 



2.2 Resummed one-loop self-energy 

While the above calculation of the interaction rate in the Born approximation is physically 
transparent, for the subsequent developments in this paper it is more convenient to obtain 
7 from the imaginary part of the self-energy. To lowest order, we can write 7(p) = 
— ImS+(p, p), with J]^{uj,p) defined as in eq. ( |1.6D in terms of the resummed one-loop 
self-energy. The corresponding diagram is displayed in Fig. ^ the blob on the photon 
line in this figure denotes the effective photon propagator of eq. (|2.13|) . 

To evaluate the one-loop diagram in Fig. ^, we use the imaginary time formalism 
and write 

m = -g'T E l^^^f^Soip-q)l.*D^'^iq). (2.17) 

In this equation, all the energy variables are purely imaginary and discrete to start with; 
namely, p^ = iun = i{2n + 1)ttT for the external fermion line, and g° = ium = i2TTmT 
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p 



p-q 



Figure 4: The resummed one- loop self-energy 



for the internal photon line, with integers n and m. Furthermore, k = p — q, So{p — q) is 
the free fermion propagator, eq. ( p..3|) , and *D'^'^{q) is the resummed photon propagator. 
We shall perform our computations in the Coulomb gauge (the one-loop damping rate is 
gauge independent [^, 0; see also Appendix B). 

The continuation of S(p) to real external energy can be done only a/ter performing 
the Matsubara sum over go = 'i^m, and consists in simply replacing (for retarded boundary 
conditions) p° = iuJn hj uj + irj, with real uj and — 0+. In order to perform the 
Matsubara sum in eq. ( |2.17] ), it is convenient to use the spectral representations of the 
various propagators. For Sq, this is given in eq. (|1.1D, with pf(p°, p) = ^Po{p^,p)- For the 
electric and magnetic photon propagators we have similarly 

'•^ dqo *Pt{qo,q) 



*At(cu,q) 
*Ai(cu,q) 



oo 27r qo — uj ' 
]_ _^ r dqo *pi{qo,q) 
(f- J-oo 2tt qo — uj 

where * pi and * pt are the corresponding spectral densities. 



*PiAqo,q) = 21m*AiAqo + ir],q) 



(2.18) 



(2.19) 



Note the subtraction performed in the spectral representation of *A;(u;,g): this is neces- 
sary since *Ai{uj,q) — 1/g^ as \uj\ oo. When the above expressions are inserted in 



eq. (|2.17|) , the sum over can be performed easily. One obtains then 

f d^q /■+°° dA;° /■+°° dg° „ , , 1 + A^(g°) - n(A;°) , 

The analytical continuation —>■ uj + ir] can now be done, and the damping rate is 
calculated as 7(p) = — ImS+(j9, p). One gets: 



7(p) 



Tig' 



UJ 



d^q /•+°° dk^ 



(2vr) 



27r J- 



do^ r 

— (5(A;° + g° - cu) 1 + iV(g°) - n{k' 
27r L 



poik) {2[uko - (p ■ q)(k • q)] * pt{q) + [cuko + (p ■ k)] *pi{q)} , (2.21 
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where u = p, = {k^, k) and k = p — q. 

The spectral functions ( |2.19| ) of the dressed photon have the following structure: 

*Ps{qo, q) = 2neiqo) z^) 5{ql - ^^(g)) + q)e{q^ - ql)^ (2.22) 

where s = I or t, Zs{q) is the residue of the time-like pole at ujs{q), and 

A(9o,g) = 37ru;p — |*A,(go,g)P, 



A(go,g) = 37cu;l^-^^^^^\*Mqo,q)\'. (2.23) 



For uj p, the energy conservation selects the positive value fco = €p-q = |p — q| from 
the spectral density pQ{ko,k) of the internal fermion. Also, the kinematics restricts the 
photon momentum to be space-like (|go| < q)- Finally, because of the infrared sensitivity 
of the damping rate, the whole contribution to F in the on-shell limit (and not only its 
divergent part) comes from soft photon momenta, q <^ T. Since, on the other hand. 



p ~ T, we can make the following kinematical approximations when evaluating eq. (|2.21| ) 
(recall that u = p): 

^P-q ^ p-p q = p-qcos9, 
uep^y - (p • q)(k ■ q) ~ p^{l - cos^ 6*), 
uko + (p ■ k) ~ 2p^, 
1 + A^(g°) - n(ep_g) - iV(g°) ~ T/g° . (2.24) 



With these simplifications, eq. ( p.21| ) becomes 



lip) - T^a'^T I I 5{qo~ qcosO) 

p,{q) + [l - cosHy p,{q)) , (2.25) 

and it is independent of the external momentum. To be consistent with the approxima- 
tions performed, we supply the above integral over q with an upper cut-off q* satisfying 
gT <^ q* <^ T. We shall verify later that, to the order of interest, the value of the integral 
is actually independent of q*. 

By using the 5-function to perform the angular integration in eq. ( ^.25| ), we obtain 

12 



In order to regularize the IR divergence, we have inserted a lower cut-off /x in the integral 
over q. Note that, because of the kinematics, the support of the energy integral is limited 
to —q<qo< q, so that only the off-shell pieces Pi^tilo, <i) of the photon spectral densities 
contribute to the damping rate. This is consistent with the physical interpretation of the 
damping rate presented in section 2.1. In fact, at this point, we can easily make contact 
between these two presentations. Namely, eq. ( p.llD in section 2.1 is essentially the same 
as the above eq. (|2.26|) , as can be seen by using eq. (|2.23| ) for the spectral densities. 
Moreover, the IR singular piece of the damping rate (|2.26|) is given by eq. ( p.l5|) , as 
we verify now through a different computation, based on the sum-rules displayed in 
Appendix A. 

Using the behaviour of these sum-rules for large photon momenta q ^ tUp, as given 
in eq. ( |A.6| ) , one can verify that 7 is independent of the arbitrary intermediate scale q* , to 
the order of interest (the contribution of the momenta g > g* is of relative order gT/q*). 
Furthermore, the infrared behavior is dominated by that term of eq. ( |2.26| ) which involves 
the transverse spectral density divided by go- Specifically, for small momenta g <^ cUp we 
can write 

r^A(go,g) = -,{l + 0{q'lul)), (2.27) 

which diverges as 1/g^ in the zero momentum limit. All the other terms give finite 
contributions as g — > (of relative order g^/cj^), and will be neglected here. By retaining 
only the leading term in eq. (|2.27|) , we obtain the singular contribution to eq. ( p.26|) : 

g^T r^^^l g^T ^^Up 



Ising = dg - = ^ In ^. (2.28) 

The upper cut-off Up ~ gT accounts approximately for the terms which have been ne- 
glected when keeping only the 1/g^ contribution to the sum-rule ( |2.27D (recall that the 
full integrand in eq. ( ^.2tj| ) is indeed cut-off at g ~ Up). As long as we are interested 
only in the coefficient of the logarithm, the precise value of this cut-off is unimportant. 
The scale Up however is uniquely determined by the physical process responsible for the 
existence of space like photons, i.e., the Landau damping. As we shall see later, this is 
the scale which fixes the long time behavior of the retarded propagator. 

In terms of collisions, the logarithmic singularity of 7, eq. ( |2.28| ), arises from the 
exchange of very soft quasistatic (go — 0) magnetic photons, as already discussed in 
section 2.1. In the present computation, this may be seen also as follows: for very soft 
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momenta q <^ Up, the function ^(go, Q')/'?o is strongly peaked at go = (see Fig. ^) and 
in the calculation of the integral ( |2.27| ) it can be replaced by the following approximate 
expression: 



1 ^ , , Svr tulq 271 



go 2 g^ + {Snuj^qo/Ay q^ 

This is, of course, just a translation of the corresponding property (|2.16|) of the magnetic 



propagator. Still, this is suggestive as it shows that the density of states Pt{%, Q')^(q'o) ~ 
(T/go)/3t(go, g) which are available for the emission (go > 0) or the absorbtion (go < 0) of 
a virtual photon with momentum q and energy go is nonvanishing in the zero-frequency 
limit go 0, in spite of the fact that the spectral density A(go,g) vanishes in the same 
limit. In fact, for very soft momenta q <ti Up, the whole density of states is concentrated 
at go = 0, as shown by eq. (|2.29| ). 



2.3 Static photon modes and non-perturbative aspects 



Because of the delta function singularity 5(go) in eq. ( p.29| ), the above discussion suggests 
that, in the imaginary time formalism, the whole IR singularity is concentrated in the 
static mode go = 0. Let us verify this explicitly by showing that, indeed, the logarithm in 
eq. ( |2.28| ) arises entirely from the magnetic contribution of the static term go = iujm = 
in the Matsubara sum of eq. ( |2.17] ) ||T^, Note that the analytic continuation of this 



term to real energy {po uj + irj) is well-defined, since all its singularities lie on the real 
axis in the complex po plane. (This is not so for the terms with go = iuJm 7^ 0, which 
individually have singularities off the real axis.) 

The magnetostatic mode gives the following contribution to the one-loop self-energy: 
S,(.;,p) = -/t|^ Y5o(^,p-q)7^*/^^^(0,q) 

2rr. [ d^g 7'V7° - (p - q) • - g'g^ sn^ 

^ J (27r)3 (cj + zr/)2 - (p - q)2 q2 " ^ ' > 

The momentum integral in eq. ( p. 301 ) shows a logarithmic ultraviolet divergence. In the 
full calculation, such a divergence would be cut-off by the contribution of the non static 
modes. (Recall the discussion after eq. ( |2.28|) .) When supplemented with an upper cut-off 



ojp, eq. ( |2.30| ) yields the following contribution to the fermion damping rate (for u; ~ p): 



i f d g i 



Ap J (27r)'^ g^ oo — p — qcosO + irj 
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~ ^ dcose6(uj-p-qcose) = aT\n , , , (2.31) 

4:71 JO J -I \UJ — P\ 

where the approximate equahty means that only regular terms have been dropped. In the 
mass-shell limit uj this reproduces the singular result of eq. ( |2.28|) . Note that the 

upper cut-off uOp is the only trace of the screening effects in the above calculation: indeed, 
the magnetostatic propagator is the same as at the tree-level, namely *A(0,g) = 

The divergence of 7 at the (resummed) one-loop level invites a closer examination 
of the higher order corrections. The two-loop self-energy is briefly discussed in Appendix 
C, where we show that the leading infrared divergence arises, again from terms where 
both the internal photons are static and magnetic. This result is readily generalized to all 
orders: the most singular contributions to the on-shell fermion self-energy are confined to 
the magnetostatic sector. When computing these contributions, all the loop integrals run 
over the three momenta q of the static internal photons, so that the infrared singularities 
are effectively those of a three-dimensional theory. Consider then a generic n-loop self- 
energy diagram with only magnetostatic modes: Its discontinuity, when evaluated on 
the tree-level mass-shell uj = p, has power-like IR divergences, possibly combined with 
logarithmic ones. Power counting shows that the leading divergences are of relative order 
{g^T / nY^^ , where /z is an IR cut-off. Such strong IR divergences are analogous to those 



identified in the analysis of the corrections to the screening mass in [0], and their presence 
signals a breakdown of perturbation theory. 

To get further insight, it is useful to consider the explicit two- loop calculation from 
Appendix C: the on-shell self-energy p) shows a linear plus logarithmic divergence. 

(There are also subleading, purely logarithmic, divergences, but these are left out in a 
leading-order calculation.) Specifically, eq. ( |C.5| ) yields 

T.t\p,p) ^ ^r(/.±(p)S(2)(p,p)) ^ .-^^ In^^ ~ ^S«(p,p), (2.32) 

where S^''(p,p) = —iaThiiuOp/jj) is the on-shell limit of the one-loop self-energy in 
eq. ( [2.31| ). Now, in order to compute the two- loop contribution to the damping rate, 
one has to expand the dispersion equation uj = p + S+(ci;,p) up to the order of interest. 
This yields the 2-loop mass-shell correction in the form 

5^(2) (p) = [z^^\p) - l]'Z^l\p,p) + (p,p) + C(3 loops) , (2.33) 
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where 



z'^^\p) - 1 



9S 



(1) 



2 aT , , 

~ (2.34) 

uj=p fJ' 



duo 

is the one-loop residue, whose leading IR-divergent part has been computed in the Ap- 
pendix. By combining eqs. (p.32|) and (|2.33|) - (|2.34|) , we note that the leading, power-like, 
divergences cancel between the two-loop self-energy and the one-loop residue, so that the 
two-loop correction to 7 is only logarithmically divergent, as at the one loop level. 

A simple argument, based on a gauge-invariant approximation to the full Dyson- 
Schwinger equation which is detailed in Appendix C, suggests that this is a general feature: 
if we assume the fermion propagator to have a simple pole at the mass-shell, then the 
damping rate remains logarithmically divergent to all orders. That is, the power-like 
divergences which occur in J]{uj = p) appear to cancel against similar divergences in 
the residue. (A similar all-order cancellation has been argued in three-dimensional QED 
at zero temperature |[45|| . ) However, the persistence of the logarithmic divergence in all 
orders of perturbation theory suggests that the analytic structure of the propagator is 
more complicated than a simple pole. 

To conclude this section, let us emphasize that when we compute the imaginary 
part of multi-loop diagrams with only static internal photons, we are actually considering 
the effects of multiple collisions involving the exchange of quasistatic magnetic photons 
with the plasma particles. The fact that these processes (or, more accurately, their most 
IR singular contributions to the interaction rate) can be effectively taken into account 
by the "dimensional reduction" to the magnetostatic photon modes is a consequence of 
the specific infrared behaviour of the resummed magnetic propagator, as expressed by 
eq. (|2l6|) . 



3 The Bloch-Nordsieck model at finite temperature 

Previously, we have shown that the leading infrared divergences in the perturbative com- 
putation of the fermion self-energy are those of an effective three-dimensional theory 
involving only static magnetic photons. We shall take advantage of this in order to get an 
explicit expression for the fermion propagator. However, before restricting ourselves to the 
static photon modes, we shall first develop a more general approach which is essentially 
a finite-temperature extension of the Bloch-Nordsieck approximation [p2| , p3 |. 

17 




Figure 6: A generic ra-loop diagram (here, n = 6) for the self-energy in quenched QED. 
3.1 Perturbation theory with soft photons 

We start by deriving a set of simphfied Feynman rules wliich allows one to compute 
the most IR singular contributions to the damping rate from higher loop self-energy 
diagrams. The leading infrared divergences arise from diagrams where all the internal 
photon lines are soft, and therefore dressed by the screening corrections. No further 
resummation of the photon lines is necessary beyond the hard termal loop approximation: 
in abelian gauge theories, all the higher order corrections to the photon polarisation tensor 
remain perturbative, and do not modify the qualitative IR behavior of the HTL-resummed 
propagator (denoted as *D^y{q)). 

Thus, when expressed in terms of the resummed photon propagator, the relevant 
self-energy diagrams contain no fermion loops: the internal photon lines are all attached 
on the incoming fermion line. A typical n-loop diagram is shown in Fig. p. There are as 
many loops as photon propagators, and we can chose all the independent loop momenta 
to be the momenta of the soft photon lines (here r = 1, n for an n-loop graph). All 
such diagrams are composed from the three following structural units: 

(i) the effective photon propagator *D^^{q); 

(ii) the fermion propagator So{p + q), where p is the hard external momentum, and q is 
a linear combination of the soft loop momenta; 

(iii) the photon-fermion vertex 7^. 

In the kinematical regime of interest, both the fermion propagator and the vertex function 
can be further simplified, along the lines explained in Section 2.2. After performing the 
Matsubara sums over the internal bosonic frequencies, and the analytic continuation to 
real external energy, the internal fermion lines are represented by spectral densities like 
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(see, e.g., eq. (^) 



Mfc°,P-q) = (A:V-(p-q)-7)^(5(fc°-ep-,)-5(fc° + 6,„,)), (3.1 



which muhiply energy denominators of the form l/{uj — kf^ — q^). Since q <^ p, we can 
use ep_g ~ep — v-q = ]9 — g cos 6* (where v = dtp/ dp = p) to replace eq. ( |3.1| ) with 

Pf (fc°, p - q) ^ (/ - p . 7)7r5(A:° - v ■ (p - q)) = /i+(p)po(A:°, p - q), (3.2) 

where the reduced spectral density 

Po(u;,p) = 27r5(cj - V • p) (3.3) 

involves only the positive-energy fermion state. The contribution of the negative-energy 
fermion state, initially present in eq. ( p.l| ), is suppressed by the corresponding large energy 
denominator. 



One sees on eq. ( |3.2| ) that neither the spin structure, nor the negative-energy fermion 
intermediate states, play an important role. In fact, the residual spin structure of eq. (|3.2| ), 
i.e. the spin matrix /i_|_(p), does not involve the loop momenta anymore, and can be 
absorbed into a redefinition of the vertex function. To see this, recall that, for a positive 
energy hard fermion, the relevant self-energy is = tr (/i_|_(p)S)/2. In the present 
kinematical regime, the spin structure of a typical n-loop contribution to S+ factorizes 
into the trace 

jmM2.../^2„ ^ itr{/i+(p)7'^i/i+(p)7^^../i+(p)7'^2n|_ (3 4) 

By using the identities (with = (1, v)) 

h^{p)Yh+{p) = v^h4p), 

tT{h4p)Y) = 2v'^, (3.5) 

one readily derives 

J^J.ll-l2 ■■■ ^J.2n _ yf^i y^J.2 y^^2n (3 6) 



The same result would have been obtained by using the reduced spectral density 
instead of ( |3.2| ), together with the effective vertex = v^. 

To conclude, the corrections to the self-energy S+ which derive from the fermion 
interactions with soft photons can be obtained from the Feynman graphs of quenched 
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QED, by evaluating the latter with the following effective Feynman rules: 

(i) the photon propagator 

(ii) the fermion (analytic) propagator 

^ . _ ^ /■+- po(fc°,p-q) ^ -1 

i-oo 27r fco-(po-go) (p^ - g") - v • (p - q) ' ^ ' 

(iii) the photon-fermion vertex F'^ = f'^. 

Any reference to the spin structure, and also to the antiparticles, has disappeared. 

We note that, when used in relation to the one-loop self-energy in Fig. 4, the above 
Feynman rules yield directly the expression ( p.26| ) for the damping rate (that is, the whole 
contribution of order g'^T, and not only its divergent piece). For higher loop diagrams 
however, we do not expect all the subleading divergences to be correctly reproduced since, 
for instance, contributions coming from mixed diagrams, where some photons are hard 
and the other ones are soft, have been ignored. 

The simplified structure which is put forward here is familiar from most treatements 



of the IR divergences at zero temperature (see, e.g., |^ and references therein). It can be 
most economically exploited within the Bloch-Nordsieck model ||3^ (see also ||33|), which, 
for the vaccuum theory, is exactly soluble. In order to search for a non-perturbative 
solution at finite temperature, we follow Ref. and reformulate this model in the 

language of path integrals. 

3.2 The Bloch-Nordsieck model in functional form 

In the Matsubara formalism, the exact fermion propagator at finite temperature can be 
obtained as 

SE{x,y) = Z-'J[dA] SE{x,y\A)exp S^-Tiln SE{x,y\A) - ^ {A,D,'a)^ (3.8) 

where SE{x,y\A) is the (imaginary-time) fermion propagator in the presence of a back- 
ground gauge field: 

-ip^SE{x,y\A) = SE{x-y), (3.9) 

and = dfj_ + igA^. In these equations, the time variables are purely imaginary (e.g., 

2^0 = -i'Tx and yo = -iTy, with < r^., Ty < (3 and 5e{x - y) = 6{t^ - Ty)6{yi - y)). 
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and the gauge fields are periodic in time, A^{t = 0) = A^(r = /?). The tree-level photon 
action has been written as 

{A,D~'A) =TY.j ^A^(^a;^,q)Z^o"i(^^™>q)^'(-^^™,-q)> (3-10) 

where Um = 27rmT with integer m, and -DoAii/ is the free photon propagator in an arbitrary 
gauge. 

The fermion propagators Se{x — y) and SE{x,y\A) are antiperiodic. For instance, 

SE{r. = 0, Ty\A) = -SE{r. = (3, r,|A), (3.11) 

and similarly for Ty. The functional determinant exp{— Tr In 5'£;(a;, y|y4)} describes the 
plasma polarization. Diagrammatically, this terms generates internal fermion loops. As 
already discussed, the only polarization effects which need to be considered are those 
contained in the photon HTL, which we denote here as 511^,^: 

TT\nSE{x,y\A) ~ ^ [a, 6UA). (3.12) 

Furthermore, the simplifications discussed in the previous subsection are easily imple- 
mented by replacing the exact propagator Se{x, y\A) in eq. ( p.8| ) with the Bloch-Nordsieck 
propagator GE{x,y\A), solution of the equation 

-t{vD,)GE{x,y\A) = 5E{x-y), (3.13) 

with antiperiodic boundary conditions. (Formally, this equation is obtained by replacing 
the Dirac matrices by the particle velocity v'^ in the full equation ( |3.9| )). With the 
above simplifications, the general equation (|3.8| ) reduces to 



SE{x,y) = J[dA]GE{x,y\A)expl^-^ [a, (3.14) 

with *-D^J = -Do/ijy + SH^iy. It is easy to verify that, when considered in perturbation 
theory, eqs. (|3.13| ) and ( |3.14| ) generate the simplified Feynman rules alluded to at the end 



of the previous subsection. 

3.3 The Bloch-Nordsieck propagator in imaginary- time 

In real-time, the equation for G{x,y\A) reads 

-t{vD.,)Gix,y\A) = 5^^\x-y), (3.15) 
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and can be solved exactly. For retarded boundary conditions, G{x,y\A) = for xq < yo, 
the solution reads 



GRix,y\A) = z^(x°-y°)5(=^)(x-y-v(x°-/))f/(a;,i/) 

dt5^^\x-y-vt)U{x,x-vt). (3.16) 



Here, U{x, y) is the parallel transporter along the straight line trajectory of velocity v 
joining x and y {y = x — vt): 

U{x, X — vt) = exp l^—ig J dsv ■ A{x — v{t — s))^ . (3-17) 



In order to verify that ( p.l6|) is indeed a solution of eq. ( p.l5|) , one may use the fact that 



the function U{x, x — vt) satisfies the following equation. 



d 

- — U{x, x-vt) = {v- D^) U{x, x-vt), (3.18) 



with the boundary condition U = 1 for t = 0. 



In imaginary-time, the resolution of eq. (p.l3 ) is complicated by the antiperiodic 
boundary conditions to be imposed on G^;: 

Gi?(r.. = 0, Ty\A) = -Ge{t, = A Ty\A), (3.19) 

and similarly for Ty. The free equation [A = 0) can be easily solved in momentum space: 

GE{lUJn,p) = (3.20) 

where ujn = (2n + 1)7tT. This is in agreement with eq. ( p.7|) . In the imaginary time 
representation, 

Ge{t,p) = Y.e-'''-^GE{^oOn,p) = e'^'P^ [^^(r)(l - n(v p)) - 0(-r)n(v p)] , (3.21) 

where n{uj) = l/(^exp(l3uj) + 1^ is the Fermi-Dirac statistical factor. 

Consider now the interacting problem, with A ^ 0. As a guidance in searching a 
solution to eq. ( |3.13| ) with antiperiodic boundary conditions, we use the solution ( p.l6| ) to 
the real-time problem, which we write in the form (with p'^ = {u, p) and v-p = uj — hfv ■ p) 

GR{x,y\A) = j -^^e-^v-i^-y)GR{xMA) 

Gr(x,p\A) = i due''^''-P^-'^'U{x,x -vt), (3.22) 
Jo 
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where the x-dependence of the function GFi{x,p\A) comes from the corresponding depen- 
dence of the background field. 

By analogy, we look for the solution GE{x,y\A) to the imaginary-time BN equation 
in the following form|^ 

GEix,y\A) = J{dp}e-'P<--y^GEix,p\A) 

GEix,p\A) = - due-''^''-P^V{x,v-p;u), (3.23) 
Jo 



where the unknown function V{x, v ■ p; m) satisfies 

V{t, = 0;u) = ViT, = P;u), 
V{x,Y-p;u = 0) + e^^^-^^V{x,Y-p;u = (3) = 1. (3.24) 

As in the real-time case, the x-dependence of Ge{x,p\A) arises entirely from its inter- 
actions with the (periodic) gauge field. If A = 0, we recover the free propagator (|3.20| ) 



by replacing V{x,v ■ p;u) with n(v ■ p), which satisfies indeed the last equation (|3.24| ) 
because of the identity 



n(e) + e^'n(e) = 1. (3.25) 



Eq. ( |3.24| ), with the indicated boundary conditions, can be solved as a series in 
powers of gA^^, that perturbative expansion: 

V ■ A{q) 
V ■ q 



V{x, V ■ p; u) = n{v ■ p) + g J [dq] e ^ ^^^'^ n(v ■ p) - n{v ■ (p + q))e 



-u{v-q) 



+ ^/MMe-^--^-)-^^:^^^:^ 

2 J V ■ qi V ■ q2 

"n(v ■ p) - n{w ■ (p + qi)) e'^^''-^^) - n(v ■ (p + qa)) e'^^''-''^) 

+n(v ■ (p + qi + qa)) e-"''-('?i+'?2)l + .... (3.26) 



It can be verified, using in particular the identity (|3.25| ), that the series (|3.26| ) satisfies 
indeed eqs. ( |3.24|) . 



^To simplify notations, the measure in the momentum integrals will be denoted below by the following 
condensed notation: 



Po,odd 
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As already noted, the quantity V{x, v ■ p; m) is the imaginary-time analogue of the 
real-time parallel transporter U{x,x — vt), eq. ( |3.17| ). This is also manifest from the 
analogy between eq. ( |3.18| ) for U{x,x — vt) and eq. ( |3.24| ) for v ■ p;m). By solving 
eq. ( |3.18| ) in perturbation theory, one generates a series analogous to ( |3.26D , where, how- 
ever, the thermal factors are absent. The correspondance between the two series can be 
easily worked out term by term. For instance. 



n{\- ■ p) — n(v ■ (p + q))e" 



-u{v-q) 



1 



and so on. In the real-time series, factorisations occur, which bring in simplifications. For 
example, in second order. 



Because of such factorisations, the real-time series corresponding to eq. ( |3.26| ) can be 
resummed into an exponential, leading to the expression ( p.l7| ). In the imaginary-time, 
the presence of the thermal factors prevents such a simple exponentiation. 



By inserting GEix,y\A), eqs. ( |3.23D and ( |3.26| ), into eq. ( |3.14| ), we can perform the 
gaussian functional integral over the photon fields term by term. This yields: 



SE{x-y) = J{dp}e-'P<^-y^SE{p) 

Se(p) = - f du e-"("-P) V'fv ■ p; 
Jo 

where V{y ■ p; u) is the functional average of V{x, v ■ p; it), eq. (|3.26| 
V(v-p;m) 

n>l 

n(v ■ p) - n{y ■ (p + qi)) e-"^^-'?^) - n{y ■ (p + qa)) e^^^^''^^) 



(3.27) 



n{y ■ P + 2^ -1 — r / dgidgs ••• dg„ — — 



+ 



... + (-l)"n(v ■ (p + qi + qa + ... + qn)) e 



-uv-{q-i+q2 + ...+qn) 



(3.28) 



and 



b{q) = ^;^*D^,(^u;„,q)^;^ 



(3.29) 



Eqs. ( p.27| )-( p.28[ ) express the Matsubara fermion propagator in the Bloch-Nordsieck 
model formal series in powers oi g . 
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3.4 The retarded propagator 



To study the mass-shell behavior of the fermion propagator, we need the retarded prop- 
agator, rather than the Matsubara one. These two propagators are related by analytic 
continuation in either the complex energy, or the complex time, plane. Here it is more 
convenient to proceed in the time representation. To this aim, we recall that the retarded 
propagator 3^(1:, p) can be obtained as 

Snit^p) = z^(t)(5>(t,p) + 5<(t,p)), (3.30) 

where the functions S-^ and S"^ are the analytic components of the time-ordered propagator E 



T3| . These can be obtained from the Matsubara propagator: 

Se{t,p) = e{T)s>iT,p)-e{-T)s<iT,p). (3.31) 

In order to get the Matsubara propagator we have to evaluate first the sum over 
Pq = iujn in eq. ( |3.27| ). Since V{y ■ p; u) is independent of po? this may be done trivially. 



by using 

T Q-i'-nir+u) ^ ^u)-5{t + u- 13). (3.32) 

n, odd 

Then, for —(3 < r < 0, we obtain 

5<(r,x) = /(^e--5<(r,p), 

^<(r,p) = e-"(^-P)V(v-p;M = -r), (3.33) 
and similarly, for < r < 

5>(r,p) = e(^-")(^-P)y(v-p;?i = /3-r). (3.34) 
In particular, the last eq. ( |3.24| ) implies 



^>(0,p) + 5<(0,p) = 1. (3.35) 

If the functions 5'^(r) and S^{t) are known explicitly, then they can be analytically 
extended in the complex time plane by simply replacing r it, with complex t. The 
functions S*^ (t) and (t) thus obtained are well defined for any t satisfying < Imt < (3, 
in the case of S^{t), and — /5 < Imt < 0, for S-^{t). For the problem at hand, these analytic 
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properties can be verified on eq. ( p.28| ): they arise from tlie fact that the thermal factors 
render the momentum integrals like 

d^g 
(2^ 

convergent for any < u < (3. We see that the statistical factors are essential to ensure 
the correct analytical properties; but, at the same time, they prevent the exponentiation 
in eq. { ^^ . 

According to eqs. ( p.30| ), ( p.33| ) and ( |3.34| ), the retarded propagator is given by 

Suit, p) = te{t)e-''^^ ■ (e'^^^ ■ P)y(v ■p;u = (3-tt)+V{v-p;u= -it)) . (3.36) 

The analytic continuation of the function V{v ■ p; u) to real time is permitted only after 
performing the Matsubara sums over the bosonic frequencies g° = ium in all the terms 
of the infinite series ( p.28|) . Fortunately, we may avoid doing this if we restrict ourselves 
to resumming the leading infrared divergences. This is further explained in the next 
subsection. 



3.5 Dimensional reduction 



In view of the discussion in section 2.3, the most IR singular terms of the perturbative 
expansion are concentrated in the static photon modes. Considering only the contribution 
of the static modes go = i^m = to eq. ( p.28| ) is equivalent to solving the Bloch-Nordsieck 
equation ( p.l3| ) in the presence of a static background field yl^(x): 



A^(x) = T 



d^g 



- e*^ ^A^i^Um 



0,q) 



T / drA^(r,x). 
Jo 



(3.37) 



With only static photon modes, the analytic continuation of eq. ( |3.28| ) to real time is 
trivial, and the sum in eq. ( |3.36| ) can be performed explicitly, term by term. As we show 
now, the thermal occupation factors compensate in this sum, and the resulting series for 
SR{t,p) can be resummed as an exponential. To be specific, consider the term of order 
(yf^ in the expansion ( p.28|) . To V^(v ■ p;u = —it), this terms contributes (-D(q) = D{0, q)) 

d^g Z)(q) 



(27r)3 (v ■ q)2 



n(vp) -n(v (p + q))e-'*(^'^) 



(3.38) 



while to e^*^^ ' P)V(v ■ p;u = (3 — it) it contributes 



d^g L)(q) 
(27r)3 (v ■ q)2 



3/3(v ■ p) 



n 



V ■ p) - n(v • (p + q)) e^(^ • q)e-'*(^ • ^^l . (3.39) 
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In the sum of these two expressions, the thermal factors disappear because of the identity 
( |3.25D , leaving 



-«^^/7l^7lM.O-— (3.40) 



(27r)3 (v ■ q)2 

By analyzing similar compensations for the higher order terms, we eventually recognize 
the power expansion of an exponential: 

Sn{t,p) = ^^(t)e~^*(^-P)A(t), (3.41) 

with 

A(t) ^ exp |-/T / (l - cost(v ■ q))} . (3.42) 

The zero-frequency photon propagator reads 

D{q) ^ *D,^{co^^ = 0, q)v-' = + 4 " ^^^] + A^^^, (3.43) 

q-^ + mjj \ q J q 

in an arbitrary gauge of the Coulomb or the covariant type (A = corresponds to both 
the Landau and the strict Coulomb gauges). The three terms in eq. ( p.43| ) corresponds 
respectively to the electric, magnetic and gauge sector. In eq. ( p.42| ), we have replaced 
the complex exponential by a cosine function, by taking into account the parity of the 
integrand. 

The g-integral in eq. (|3.42|) presents a spurious ultraviolet logarithmic divergence in 
the physical sector (i.e., for electric and magnetic photons). This divergence is unphys- 
ical since in the full theory, including also the non-static photon modes, the g-integral 
would be cut-off at momenta g ~ cUp (recall the discussion in section 2.2). Thus, to be 
consistent with the approximations performed, we have to complement the above "dimen- 
sional reduction" with the prescription that an upper cut-off of the order gT is added in 
momentum integrals, in the physical sector. Since this cut-off is not exactly known, it 
will be important in what follows to verify that the physical predictions are independent 
from its precise value. In the gauge sector, on the other hand, no such cut-off is needed 
since the corresponding momentum integral turns out to be ultraviolet finite (see below, 
eq. (CT) ). 



Eq. ( |3.42|) determines the large time behavior of the fermion propagator, to be dis- 



cussed in the next section. At a first sight, the considerable simplifications leading to 
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this equation (and coming from the restriction to the static photon modes in eq. ( p.28| )) 
may seem rather accidentaL However, as we explain now, there is a simple reason for 
these simplifications, and, in fact, eq. ( p.42| ) could have been obtained in a more direct 
waypi[|, which avoids some of the complications of the Matsubara formalism (the latter 
are essential only for the non-static modes). Let us indeed return briefly to the functional 
integral of eq. (|3.14| ), and consider its approximation where the Bloch-Nordsieck propaga- 
tor Ge{x, y\A) includes only the static electromagnetic field A^(x) of eq. ( p.37|) . Then, the 
contribution of the non-static photon modes to the functional integral trivially factorizes, 
and is compensated by the corresponding contribution to the partition function Z, thus 
leaving 

SEix, y) = Z,' I [dA] GEix, y\A) exp {a, *D-'a)J^ , (3.44) 



where = Afj_{uJm = 0, q), and {A, *D~^Aj^ denotes the Um = contribution to the effec- 
tive photon action ( |3.1CI| ); correspondingly, Zq is the partition function of the static mode 
alone. Since the background field ( |3.37D is time- independent, the propagator GE{x,y\A) 
depends only on the time difference Xq — yo, i-e. GE{x,y\A) = Ge{xq — yQ,^,y\A). Its 
Fourier transform can be analytically continued in the complex energy plane, and the 
resulting function coincides, in the upper half plane, with the retarded propagator. It 
is then convenient to take the Fourier transform of eq. ( p.44| ), and write [p^ = {iu)n,p), 
Un = (2n + l)7rT): 

Se{p) ^ Z,' J [dA] Ge{^,p\A) exp |-i [A, J . (3.45) 

Since the energy pq enters eq. ( p.45| ) as an external parameter, the continuation to real 
external energy po uj+ir], and the Fourier transform to real time, can both be performed 
before doing the functional integration. Thus, the retarded propagator Sji{t,x.) can be 
directly obtained as the functional average of Gr{x, y\A), which is known explicitly (recall 
eqs. ( gl6D and (|3l^)). 

Specifically, eqs. (|3.45|) and ( p. 161) give S'ij(t, p) in the form (|3.41|) , where 

A(t) = Z,' J [dA] U{x, X - vt) exp |-i (A, *D-'a)^^ , (3.46) 
and the parallel transporter is that of a static background field: 

U{x,x-vt) = exp|-yd=^?/j^(y)A^(y)|, 

jM = igvj'ds6^'\^-y-vs). (3.47) 

•/ 
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A straightforward computation yields then 

A(t) = exp|iry'd3xid3x2j'^(xi)*D^,(xi-X2)/(x2)} 



exp 



T / dsi / ds2V{^r{s^-S2)) 



(3.48) 



where 



r>(x) 



(3.49) 



By using the last equation 



is the Fourier transform of the static photon propagator 
to perform the Si and S2 integrations, we may cast eq. (|3.48|) in the form (|3.42|). 



4 The infrared structure of the fermion propagator 
4.1 Large time behavior 

The non-trivial time dependence of the fermion propagator is contained in the function 
A(t), eq. ( |3.42| ). Because our approximations preserve only the leading infrared behavior 
of the perturbation theory, eq. ( p.42[ ) describes only the leading large-time behavior of 



A(t). Since the only energy scale in the momentum integral of eq. ( |3.42| ) is the upper 
cut-off, of order gT, the large-time regime is achieved for t ^ l/gT. 

The expansion of eq. (|3.42[) in powers of g^ reproduces the dominant singularities of 



the usual perturbative expansion for the self-energy. Let us verify this for the correction 
of order g"^: 

SSd^, p) = -/r , /" dt e«<"-v • P+".) / (1 - COS i(v . q)) . (4.1) 

We perform first the time integration and obtain, after simple algebraic manipulationsj^ 

6Sr{uj, p) = -Go{iu, p)S(u;, p) Go{iu, p), (4.2) 

where 

Go(uj,p) = I / dte'*("-^ = -, (4.3) 

JO — V ■ p + «?7 

^ The self-energy S which appears here corresponds to the spin projection S+ of the fuU self-energy. 



(See eq. (jlj).) 
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is the free BN propagator, and 

E(.;,p) = -/r/^D(q) ~\ (4.4) 

J [2-nY — V ■ (p + q) + zr] 

The imaginary part of this equation determines the damping rate according to 7 = 
— Im Tj{ijJ = p). We can write, with e = — v ■ p. 



ImS(^,p) = -Trg^T J -—-^6{u -V ■ {p + q))D{q 



i2ny 

^^"^^^d.Ji,fl-£;i-^.A^|, (4.5) 



Att J\e\ \ t) q +mj-) q 

which, in the mass-shell limit e — 0, and with an IR cut-off /x in the magnetic sector, 
yields (with a = g^/4:TT) 

This first piece inside the parantheses, which comes from the magnetic sector, reproduces 



the singular piece of the resummed one- loop calculation (recall eq. ( p.28|) ). On the other 
hand, the other two pieces are not correctly reproduced by the present calculation. The 
electric piece, which is finite and of the order g^T, occurs even with a minus sign (recall 



that the contribution of the electric scattering to the interaction rate in eq. (|2.26| ) was 



positive). The gauge- dependent piece turns out to be non- vanishing, but it could be 
eliminated by introducing an IR cut-off fi in the gauge sector as well, and by taking 
the on-shell limit only subsequently (see also the discussion in Appendix B). This 
situation is generic: our approximation yields correctly only the leading IR divergences 
of the usual perturbation theory, which all arise from the magnetic sector, but not the 
subleading terms. In particular, the contributions involving the electric and the gauge 
sector are subleading, and should be discarded for consistency. This is equivalent to using 
D{q) = 1/q^ rather than the full static propagator of eq. (|3.43| ). 



Let us verify now that the full, non-perturbative, expression of A(t), eq. (|3.42|) , 
is free of infrared singularities. Inspection of the integrand in eq. ( p.42|) shows that the 
dominant IR behaviour arises from the limit |v ■ q| = qcosO 0. This is consistent with 
the calculations in sections 2.1 - 2.2 showing that the divergences come from the exchange 
of magnetic photons emitted or absorbed at nearly 90 degrees. We have, in this limit. 



2 



1 — cost(v ■ q) t 

(T^lip ~ 2"^ 
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0(t^(vq)2), (4.7) 



and the momentum integral is IR safe, as advertised. We see here, once again, that 
the gauge dependent piece of the photon propagator ( |3.43| ) does not contribute to the 
leading IR behavior (which is given by the term in of the magnetic propagator). 
Indeed, because of the factor (cos6')^, the gauge propagator {cosOY /q^ is less singular as 
qcosO 0. 

Consider now the UV behavior of the g-integral. This depends logarithmically on 
the UV cut-off ~ Up, and, as a consequence, the large time behavior of A{t) is insensitive 
to both the precise value of the UV cut-off, and to the specific procedure which is used 
for its implementation. This will be verified explicitly below. 

The evaluation of A(t) is most simply done by using the coordinate space represen- 
tation ( p. 4^ ) for A(t). Corresponding to D{q) = we have T>{x.) = l/Anx, and we 
obtain, for t > l/up, A(t) = exp{-g^T F(t)), with 

F{t) = i ^*dsi^*ds2:P(v(si-S2)) 

= ^Vd./rd..^l^if^^^^.f (ln.,t+0(l)). (4.8) 
Sir Jo Jo pi — S2I 47r 

In this calculation, the ultraviolet cut-off has been introduced in the function 6'(|si — S2I — 
l/wp). Let us verify that the same large time behavior is obtained with a different UV 
regularisation, namely, with the modified photon propagator D{q) = 1/q^ — l/(g^ + uj^) 
(Pauli-Villars regularisation). By using 

we get successively 

F{t) = — dsJ dS2- 



57r Jo Jo \Si — S2I 

t f 1 - e"""" * r^f * , 1 



1 + ds 



Air [ ujpt Jo ^ ) 

= f |ln^,t + (^^ - 1) + ^ - ^M-^P^) + E^iupt)] , (4.10) 

4:11 [ UJpt J 

where Ei(a;) is the exponential- integral function, Ei(x) = dy {e^^"^ / y) , and •je the 
Euler constant. At very large times, Upt ^ 1, we may use the asymptotic expansion of 
the exponential- integral to get, for the r.h.s. of eq. ( [4.10|) , 

F{t) ~ {t/An)(lnu;pt+i-fE - 1)) ^ (t/47r) IncUpt, (4.11) 
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which coincides, as long as the leading logarithm is concerned, with the previous result 
( [4.8[ ). On the other hand, the subleading term, i.e. the constant under the logarithm, is 
dependent on the UV regularisation. Thus, as expected, it is only the dominant behavior 
at very large times which is consistently described by our approximation; the subleading 
terms should be ignored. 

We have argued before that the gauge-fixing terms are not important to the order of 
interest. To verify this explicitly, we compute the gauge-dependent contribution to F{t), 
as given by the last term of the photon propagator (|3.43|) : 

SFit) ^xfA-. ^ '"r^^^T^^ = A ^. (4.12) 
^ ^ J (27r)3 g4 ■ q)2 Svr ^ ^ 

At large times, this is indeed subleading with respect to ( 4.11 ). Note that, the momentum 
integral in eq. ( [4.12|) being ultraviolet finite, no upper cut-off has been necessary in its 
evaluation. 

We conclude that, at times t ^ 1/t^p, the function A(t) is gauge-independent and 
of the form {a = g'^/Air) 

A{ujpt > 1) ~ exp(-aTtlna;pt). (4.13) 

The most striking feature of this result is the fact that, at very large times {upt oo), 
the fermion propagator is decreasing faster than any exponential. We also note that the 
scale of the time variations is fixed by the plasma frequency Up ~ gT. 

A measure of the decay time r is given by 

- = «Tlnu;„r = aT('ln^-lnln^+ ..) . (4.14) 
r V aT aT J 

Since aT ~ gUp, we see that r ~ 1 / {g'^T \ia.{l / g)) . This is very close to the perturbative 

result ( p. 281 ), which, in the presence of an IR cut-off ~ g'^T, predicts a damping rate 

7 ~ g'^Tln{l/g). A comparison of the two decay laws, A(t) = exp{—g'^TF{t)), with F{t) 

from eq. ( |4.1CI| ), and the exponential|U] Ai;,(t) = exp(— 7t) with 7 = aTln(l/5f), is presented 

in Fig. 1^ for (jf = 0.4. In this figure, the time is measured in units of l/ujp, and the results 

displayed for A(t) can be trusted for values Upt » 1, where our approximations are 

expected to hold. For very large times, t ^ t, the function A{t) is indeed more rapidly 

decreasing than the exponential Ailt). On the other hand, for intermediate, but still 

II This is the spectral function which would produce an exponential decay in time with a lifetime as 
close as possible to the nonperturbative result (4.14). 
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Figure 7: The time behavior of the fermion propagator as described by the nonperturba- 
tive result A(t) (full line) and by the exponential Ai(t) (dashed line) for g = 0.4. On the 
abscissa axis, time is measured in units of 1/ujp. 
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large, times, 1/gT ^ t <^ l/g^T, the opposite situation holds: A(t) > /S.L{t). When 
discussing the lifetime of the excitation, it is rather the intermediate range of times which 
matters, since for asymptotically large times t ^ 1/ g'^T the excitation has already decayed. 
It follows that, for the range of times of interest, the decay of the excitation is actually 
slower than the one predicted by perturbation theory. 



4.2 Mass-shell behavior 

The non-trivial large-time behavior exhibited in eq. ([4.13|) has interesting consequences 
on the behavior of the retarded propagator in the complex energy plane. In fact, since 
at large times A(t) is decreasing faster than any exponential, the time-integral giving the 
Fourier transform 

/oo roo 
dte-^-%(t,p) = W dte^*(--^P+^'')A(t), (4.15) 
-oo Jo 

is absolutely convergent for any complex (and finite) u. That is, the retarded propagator 
Sii{uj) is an entire function, with sole singularity at Imu — oo. Recall, however, that 
strictly speaking, our present approximation holds only in the vicinity of the mass-shell. 
Therefore, when speaking about — v ■ p| — > oo we have in mind off-shell deviations 
which are much larger than g^T. To further clarify this point, let us give a crude estimate 
of how Sr{u!) increases as Imtu — — oo. To this aim, let us consider a; = v ■ p — i^, with 
real and positive (. We write: 

POO 

SR{0=t dte^'A{t), (4.16) 







which is a purely imaginary function of (, and consider the behaviour of |S'ij(C)| for 
( ^ aT. Regarded as a function of t, the integrand e''*A(t) is rapidly increasing for 
small t, but it is decreasing for sufficiently large values of t, where the decay of A(t) starts 



to dominate. Assuming the time integral in ([4.1 6|) to be dominated by large values of t, 
— which is correct for large enough ( — , we can use the asymptotic expression ( |4.13|) , 
and determine the time t* at which the integrand is maximum: 

ujpt* = exp^^^. (4.17) 
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By using the fact that the integrand is positive definite, and that, according to eq. ( |4.17| ), 
C — aT In ujpt* = aT and thus ( — aT In Upt > aT for any t <t*, we can writeQ 

ri* 1 



SO that 



iSniOl > dte"^* = — (exp(arr)-l), (4.18) 



1^^(01 > ^ jexp (^exp^) - lU ^exp (^exp^) , (4.19) 



where g = aT/up = {S/Air) g. Eq. ( [4.19|) shows that |S'ij(C)| is rapidly increasing starting 
with values of ( such that ge'^ ~ 1, that is, ( ~ aTln(l/5f). In perturbation theory, Sr^u) 
has a pole at = v ■ p — ^7, where 7 ~ aTln{ujp/ fi) ^ aTln{l/g) if yU ~ g'^T. Thus, our 
non-perturbative solution for Sr^u) replaces the pole at finite distance by an essential 
singularity at —zoo, which however starts manifesting itself at distances ~ g'^Tln{l/g) 
below the real axis, that is, at the same distances as the pole of the perturbation theory. 

Since Sr{uj) is analytic in any finite neighbourhood of the tree-level mass-shell at 
cij = V ■ p, we need to clarify the mass-shell interpretation. To this aim, we consider the 
spectral density p{u, p) 

POO 

p{uj, p) = 2 Im Sr{uj, p) = 2 / dt cos t{v ■ p) A(t), (4.20) 

Jo 

where v ■ p = u — v ■ p. It satisfies the sum-ruleQ 

du 



r-^p(a;,p)=A(t = 0) = l. (4.21) 

J-oo ZiT 



We have calculated p{uj, p) numerically, and the result is plotted, for a coupling constant 
g = 0.08, in Fig. p. We also represent, for the same value of g, the Lorentzian spectral 
function [e = v ■ p) 

PL(e) = ^5^, (4.22) 
+ 7^ 

with 7 = aTln{l/g). This is the spectral function which would produce the exponential 
time decay ^L{t) = exp(— 7^) alluded to at the end of the previous subsection. It is 
seen on these figures that, in the weak coupling limit, the spectral density p(e) has the 
**This estimate was suggested to us by A. Rebhan [|2| 

''"I'ln fact, this sum rule holds exactly in the Bloch-Nordsieck model, independently of the restriction to 
the static photon mode. In general, A{t = 0) is replaced, in eq. (4.21), by S^{t = 0, p) + S'^{t — 0,p), 
which is also equal to one, as shown by eq. (3.35). 
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Figure 8: The spectral density p(e) (full line, in units of l/ujp) and the lorentzian pL(e) 
(dashed line) for g = 0.08, as a function of e = f ■ p in units of ujp. 

shape of a resonance strongly peaked around e = 0, and with a typical witdh of the order 
1/r ~ g'^Tln{l/g), that is, of the same order as that of the Lorentzian. This allows us 
to identify the mass-shell of the full propagator at = v ■ p, as at treel-level. Moreover, 
it is clear from Fig. ^ that, for very small g 1, the nonperturbative spectral density 
is even sharper than a Lorentzian. Thus the net result of the infrared effects considered 
here is to slightly enhance the stability of the quasiparticle state (see also the discussion 
at the end of the previous subsection). 

Finally, it is interesting to compute the imaginary part of the exact self-energy, by 
inverting the Dyson- Schwinger equation S]i^{tu,p) = —{u — v • p) + J2r{uj,p). A simple 
calculation yields 

where e = v ■ p, p(e) is the spectral density of eq. ( [4.20|) and 

/•oo 

(T(e) = 2 / dt sinet A(t). (4.24) 
Jo 

This is represented graphically in Fig. ^ together with the pure one- loop result, eq. ( [4.5| ), 
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Figure 9: The imaginary part of the self-energy, eq. (|4.23| ), as a function of the energy, 
for g = 0.08: nonperturbative calculation (full line), one-loop result (dotted line) and 
one-loop result in the presence of an IR cut-off ~ g^T (dashed line). All the quantities 
are measured in units of ujp. 
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which shows a logarithmic divergence as e — (dotted hne), and the screened one-loop 
result, as obtained from eq. ( [4.4[ ) after inserting an IR cut-off equal to aT (dashed line). 
As manifest on this figure, the full result for Im is finite at the mass-shell e = 0, and 
inferior to the value predicted by the perturbation theory with an IR cut-off ~ g^T. The 
latter property is consistent with the the previous analysis of the spectral density, and 
also of the time behavior at intermediate times. One can also verify the non perturbative 
character of the solution. For example, ImE/j(e = 0) = — 1/ dt A{t) has no expansion 
in powers of g"^ even if one keeps ujp constant in eq. ( [4.13| ) for A(t). 



5 The lifetime of the soft fermionic excitations 



For soft momenta, p ~ gT, the quasiparticles become collective excitations, with non- 
trivial dispersion relations Q and self- interactions [^, To leading order in g, the 
dispersion relations are real, and the quasiparticles propagate without damping. At next 
to leading order, coUisional damping occurs. The corresponding damping rate 7 has 
been calculated in the effective (i.e., HTL-resummed) perturbation theory]^. For an 
excitation with zero momentum (p = 0), 7 is finite and of the order g'^T 0, ^f. However, 
for excitations with finite momentum p ^ g^T, the lowest order perturbative calculation 
of 7 meets with the same infrared problem as that discussed for the hard particles [|18|, [25 



As we shall see, this problem is solved by the same technique as that used for the hard 
fermion. 



5.1 The HTL approximation 

Let us recall first the main features of the dispersion relations for soft fermions, to leading 
order in g. They are obtained from the poles of the effective propagator *S{u!, p) which 
is obtained as *S^^ = Sq'^ + with ST,(uj, p) denoting the fermion self-energy in the 
HTL approximation Jl], |2|: 

SE{co,p) = coUf' ^— -• (5.1) 

J Att 00 — v ■ p + trj 

In this equation, uiq = gT/^/S is the frequency of the spatially uniform (p = 0) fermionic 
excitations. The propagator is conveniently written in the form ( |1.7| ), that is, 

*S{u;,p) = *A+(a;,p)/i+(p) + *A_(^,p)/i_(p), (5.2) 



38 



where 

*A^{u;,p) = ~^ --, (5.3) 

LOT (p + ()S±(a;,p)) 

and 

5S±(a;,p) = ±itr(/i±(p)5S(^,p)). (5.4) 

The pole equations *A'^^{uj{p),p) = yield two positive energy branches uj±{p) instead 
of the usual one (with u = p) in the free electron spectrum. For u close to the mass-shell 

at ujs{p), s = ±, we can write 

*As(^,p) ^ "-f^ , (5.5) 

where Zs{p) is the residue of the mode s , 

, _ d6Es{u,p) 

Since uj±{p) > p for any p, both dispersion relations are real: the quasiparticles propagate 
without damping in this approximation. For small momenta, p <^ gT, uj±{p) ~ ujq ±p/3. 
The upper branch is strictly increasing (f+(p) > for any p), while the lower branch has 
a minimum at p = ~ 0.92 uq. At very large momenta, p ^ ujq, both branches approach 
the light cone, but z+(p) — * 1, while z_{p) vanish exponentially. (See Refs. |jlO|, [l^, |l^ 
for more details and physical interpretation.) 

Because of the gauge symmetry, the nonlocal character of the HTL self-energy (|5.1| ) 
leads to effective interactions between a fermion pair and any number of soft photons. 
For instance, the Ward identity 

*T,{p,p + q) = *S-\p) - *S-\p + q) , (5.7) 

requires the existence of a nonlocal 3-point vertex function, which is indeed found in the 
form T^(p,p + g) = + 5r^(p,P + <?), where 5T^j_{p,p + q) is the 3-point HTL [|, |]: 

/liP^P + l) j Att {v ■ p + ir]){v ■ {p + q) + irf) ^ ^ 

Similarly, higher vertices, without analogue at the tree-level, are necessary in order to 
fulfill the higher order Ward identities. We show here one more example, namely the 
Ward identity satisfied by the 2-fermions — 2-photons vertex function: 

9i *r^;.(pi,P2;gi,g2) = *^v{.PuPi + ^2) - *^u{.pi + + qi + ^2) , (5.9) 
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(a) (b) 

Figure 10: One-loop diagrams for the soft fermion self-energy in the effective expansion. 

where, in the left hand side, pi and qi are respectively the momenta of the incoming 
fermions and photons, with pi + P2 + Qi + (I2 = 0. For what follows, it is important to 
remark that all the HTL vertex functions are (almost) uniquely determined by the self- 
energy ( |5.1| ) and the Ward identities like eq. ( |5.7| ) 0, This is so since the non-linear 
structure of the effective action of the HTL's is the minimal one which is consistent with 
the gauge symmetry]^, |TT[ . 



5.2 Perturbation theory for the damping rate 

In this section, we discuss the perturbative computation of the damping rate for the 
soft fermion, and the related IR problems. After a brief summary of the leading-order 
computation 0, 18, 25 1, we discuss higher orders corrections and how they simplify in the 



computation of the leading divergent terms. 

The dominant contribution to the damping rate, of order g^T, comes from the 
imaginary part of the (resummed) one-loop self-energy *S(co',p), as given by the two 



diagrams in Fig. |T0| [Q. Specifically, 

7±(p) = -z±{p)lm*T,±{uj±{p) +iT],p), (5.10) 

where *T,±{p) = tr (/i-j-(p)*S(p))/2 and the subscripts ± refer to the two positive-energy 
modes in the fermion spectrum. Note that, in general, the "tadpole" diagram in Fig. |T^.b 
gives a non-trivial contribution to 7, since the 4-point vertex itself has a non-zero discon- 
tinuity. Moreover, the imaginary part of the diagram in Fig. |r^.a comes not only from 
the cutting of the internal propagators (as for the usual one-loop diagram discussed in 
section 2), but also from the discontinuity of the resummed 3-point vertices. 

In what follows, we concentrate on the singular contribution to 7. This comes 
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entirely from the diagram in Fig. ITO.a |T8|, which reads 



^S,(p) = -^?2t Y: I S\*^MP + <iyS{p + qrT,{p + q,prD^''{q). (5.11) 



It has been aheady recognized PB| that the singular piece of 7 arises from the same 
kinematical regime as for a hard fermion, namely, from the exchange of a very soft (g ^ 
g'^T) magnetic photon at nearly 90 degrees. This allows for kinematical approximations 
identical to those encountered in section 2. In particular, the whole singularity can be 



reproduced by restricting the calculation to the magneto static mode q = icOm = p5 
with propagator *D^^{uJm = 0, q) = 6^^ /q"^, and with an upper cut-off ~ gT. Furthermore, 
the internal fermion propagator *S'(u;,p + q) is nearly on-shell, since u ~ i^±{p), and 
q <^p. Thus we can write co'±(p + q) ~ t^±(p)+v±(p)-q, — where \s{p) denotes the group 
velocity of the mode s, Vs(p) = duJs{p)/dY> = Vs{p)'p — , and replace *5'(co',p + q) — > 
*A±(u;, p + q) /i±(p), with (recall eq. (|5.5|) ) 



*A±(^,p + q) ^ ^^^^^ , (5.12) 

uj - uj±{p) - v±(p) ■ q 

where the upper (lower) sign applies according to whether the external line is close to the 
mass-shell of the upper branch, or of the lower branch, respectively. A final simplifica- 
tion refers to the 3-point vertex function *r*(p, p + q), where we can neglect the photon 
momentum q and use the differential form of the Ward identity ( |5.7| ) to write: 

•r-(p.p) ^ (5.13) 

(The dependence of the vertex function on the external energy uj is not indicated explic- 
itly.) The inverse propagator is conveniently written as (recall eq. (|1.7| )) 

*S-\lu,p) = *A-\uj,p)h4p)+*Az\iu,p)h4p). (5.14) 

From eqs. ( p.l3| ) and (pTTj), we obtain, for u ^ ^^±{p), 

ft±(p)T'(p,p)ft±(p) = ft^(p) ^-i^ A^(p) -±|W 

tr(/,.(prr'(p.p)) = 2^^2|M. (5.15) 

The particular spin projections of *r' written down above are the only ones which enter 
= tr {h±. *S) /2, and therefore the damping rate (|5.1CI|) . 
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Note that the simphfied vertex ( |5.15| ) has no discontinuity, so that the whole imagi- 
nary part of *S in the kinematical regime of interest arises by cutting the internal lines in 
Fig. [ly.a. Specifically, the previous approximations yield the dominant (infrared singular) 
piece of the one-loop damping rate as[]TB|, ^| 



7±(P) - z± g^T [ -^-^ ^ ^ — Im 

J 27r z-i- z+ oj — 



-z± 



(27r)'^ g2 uj — uj±[p) — \v±\q cos 6 + ir] 



aT\v±{p)\\Ti- (5.16) 



which is very close to eq. (|2.31| ) for a hard fermion (recall that |v| = 1 for the hard 
quasiparticle). 

Consider now the higher order corrections to 7, with emphasis on the leading in- 
frared contributions. By relying mostly on the gauge symmetry, we argue now that the 
most singular contributions to 7 arise from multi-loop diagrams which involve the (re- 
summed) 3-point photon-fermion vertex, but not the higher order vertice^. This is so 
since in the kinematical regime of interest, the inverse fermion propagator, 

*A;^(^,p + u,{p) - v,(p) ■ q) , (5.17) 

is linear in the photon momentum q, so that the Ward identity (|5.7|) can be satisfied by 
a 3-point vertex *r*(p, p + q) which is independent of the momentum of the photon leg. 
And we have seen indeed that the singular one-loop contribution is obtained by replacing 
*r*(p, p + q) with *r*(p, p), which is independent of q and (up to a spin projector) equal 
to vKp)/ Zs{p). Furthermore, with a q-independent 3-point vertex, all the other, higher. 
Ward identities — as the one shown in eq. ( |5.9| ) — are trivially satisfied by setting the 
n-point HTL's with n > 4 to zero. Since, as alluded to before, the vertex HTL's are 
essentially determined by the Ward identities, it follows that the higher-point vertices 
(beyond the 3-point function) are not important in the kinematical regime of interest. 

We thus conclude that, in order to isolate the most singular contributions to 7^ 
[s = ±) in perturbation theory, we have to consider the same diagrams as for the hard 
fermion, and evaluate them with the following simplified Feynman rules: 

(i) the photon propagator DQ^q) = 5''^ /q^; 

(ii) the fermion propagator */S.s{uj, p + q) from eq. ( ^.12|) ; 

"This can be also verified by power counting, as in Ref. for the one-loop calculation. Namely, 
cutting a vertex rather than a fermion propagator, yields a factor of l/(v<;(p) • q) less, and thus a less 
singular infrared behavior. 
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(iii) the photon-fermion vertex *r*(p) = vl{p)/zs{p). 

The momentum integrals over the photon momenta should be computed with an upper 
cut-off of the order Up. Strictly speaking, the above simplifications hold only for very soft 
momenta, q <^ p ^ gT, and not up to momenta q ~ oOp. This is not important, however, 
since the dominant (singular) contributions arise from the limit g ^ and are insensitive 
to the upper cut-off. 

Note that the above Feynman rules are essentially those of a local effective field 
theory, in contrast with the general HTL Feynman rules, which are non-local. (The 
apparent dependence on u and p is irrelevant here, since these are the fixed energy and 
momentum of the external line; they enter the computation as parameters). Furthermore, 
the IR contribution to 7 is largerly insensitive to the details of the HTL resummation, 
which enters only via the global factors Vs{p) and Zs{p). Actually, to the order of interest, 
7s is even independent of the residue Zs{p), as also suggested by the one- loop result ( |5.16|) . 
This is so since a general n-loop graph contributing to (in the simplified perturbation 
theory introduced above) involves 2n vertices *r*, and therefore a factor z~'^^, and (2n — 1) 
propagators *As, which yield a factor zf^'^. The remaining factor of 1/zs disappears in 
the computation of 7^ = —ZshxvTjs. 

5.3 The Bloch-Nordsieck model for a soft fermion 

At this point, the analysis of the dominant mass-shell behavior of the soft fermion becomes 
almost identical to the corresponding analysis for the hard fermion. This this analogy 
is due to the fact that the soft photons responsible for the IR divergences have typical 
momenta q <^ gT, which are much smaller than the momentum p ~ gT of the soft 
fermion. In view of this, the whole discussion in sections 3 and 4 can be directly extended 
to the case of a soft fermion. 

Specifically, the simplified Feynman rules which apply in the IR regime are, once 
again, those of the Bloch-Nordsieck model, and can be summarized in the following func- 
tional integral representation of the soft fermion propagator: 

5±(x, y) = Z,' I [dA] G^ix, y\A) exp |-i (A, D,'A)Jj . (5.18) 

In this equation, Gs{x,y\A) is the Bloch-Nordsieck propagator for the quasiparticle in the 
mode s, s = ±, in the presence of the static magnetic field = (0, A(x)), and satisfies 
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(with < = (l,v,)) 

-i{vs-D,)Gs{x,y\A) = Zs5^'^\x-y). (5.19) 

Furthermore, 

(A,Do-ia)^ = 1/ d'xd?yA\^)D,y{^-y)A^{y), (5.20) 



where the vector field y4*(x) has been defined in eq. ( p.37|) , and Doij(x.) = Sij/Anx. Note 



that the free (retarded) BN propagator, as obtained from eq. (|5.19|) with A = 0, reads 



G^iu, p + q) = \ . (5.21) 

u; - Vs • (p + q) + z?7 

Strictly speaking, the mass-shell for the BN particle of momentum p, that is = Vg ■ p, is 
different from the real leading-order mass-shell, at u = ujs{p). This is so, of course, since 
the dispersion relations for soft fermions are not linear, so that the group velocity |vs| 
is really momentum dependent. However, this difference is not important, since the BN 
propagator ( |5.21| ) presents the correct dependence on q in the mass-shell limit. Compare 
in this respect eqs. ( [5.21[ ) and ( p.l2[ ): in both these equations, it is the difference in 



energy with respect to the mass-shell which matters, rather than the precise value of the 
mass-shell energy itself. For u = ujs{p) in eq. (|5.12|) , and respectively for = Vg ■ p in 
eq. (|5.21|) , the propagators in these two equations become identical. 



Eqs. ( |5.18|) -( ^.19|) are further manipulated as in section 3.5 (recall, especially, the 
discussion after eq. ( p.44| )). As a result, we obtain the retarded propagator for the two 



fermionic modes ±, for momenta p ~ gT and energies close to the mass-shell, uj ^ uj±{p). 
It reads 

roc 

S±iuj,p) = tz^ip) dte^*(--'^±(^)+^^)A±(t), 

Jo 

A,(t) = exp |-,^T / ^^^^ (l - cost(v, ■ q))} . (5.22) 

In this equation, D±{q) = t>!|_Z}Q (q)t;j_ = f±/g^, so that we can write 

A±(t)= A{\v±\t), (5.23) 

with A{t) as given by eq. ( p.42| ) where -D(q) — > 1/g^ and v is an arbitrary unit vector. 
Note that the functions A± (t) are implicitly dependent on the momentum p, via the group 



velocities v±{p). Both the mass-shell behavior of the propagator (|5.22|) and the large time 
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behavior of the propagator S±{t,p) follows from the analysis in section 4. At very large 
times LJp\v±\t ^ 1, we have 

A±{ujp\vj-\t 3> 1) ~ exp [-aT\v±\t\n{iUp\v±\t)y (5.24) 

The spectral density of the mode s is peaked around uj = ujs{p), with a width of the 
order g'^T\vs\lTi{l/g). In particular, for the lower mode uj^{p), and for p = pc, where 
V-{pc) = 0, eq. (|5.23| ) shows that, to this approximation, the "plasmino" mode is not 
damped, in accordance with the one- loop result for the damping rate, eq. ( ^.16| ). 

6 Conclusions 

The analysis presented in this paper suggests that the damping of the fermionic excita- 
tions with momenta p ^ g^T is not exponential in time, but of the more complicated 
form Sn{t) ~ e"'^'-^-'* exp{— a;T|f 1 1 ln(co'p|f |t)}, where v = dE/dp is the group velocity 
of the excitation, Up ~ gT is the plasma frequency, and a = (7^/47r. As a consequence, 
the retarded propagator Sr{uj) has no quasiparticle pole, but the spectral density shows 
nevertheless a sharp resonance peaked at = E{p), with a width ~ g'^T ln{l / g) . At the 
present level of accuracy, the mean energy E{p) is given by the leading-order approxima- 
tion, namely E{p) = p for a hard excitation, and E{p) = ujj-{p) for a soft one. We note that 
this result solves the IR problem of the damping rate in a very "soft" way, by essentially 
replacing the IR cut-off fi in the perturbative result Ai(t) = exp{— aT|t>| t \n{ujp\v\/ fi)} 
with the inverse of the time. Thus, quantitatively, the lifetime of the excitation does not 
differ much from that obtained from leading-order perturbation theory. 

The asymptotic behaviour of the retarded propagator has been obtained by solv- 
ing exactly an effective theory which reproduces all the leading infrared divergences of 
the perturbation theory. The physical processes which are responsible for these diver- 
gences are the multiple collisions involving the exchange of long wavelength, quasi-static, 
magnetic photons, which are not screened by plasma effects. By comparison, the longitu- 
dinal, gauge sector is less singular in perturbation theory, and does not contribute to the 
dominant large time behaviour of the non-perturbative solution. 

At finite temperature, the presence of the thermal bath amplifies the IR divergences, 
in such a way that they become effectively those of a three-dimensional gauge theory. 
Then, a comparison with massive QEDg |^ helps explaining why an IR divergence occurs 
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for the one- loop damping rate, in contrast to the zero temperature case where the IR 



problem does not affect the dispersion equation, but only the residue of the propagator [Q. 
At this point, we should recall that the explicit solution that we have proposed here relies 
essentially on the 3-dimensional character of the dominant singularities. This has been 
widely recognized in relation with the infrared structure of thermal field theories |^ , and, 
in the calculation of static quantities (like the free energy or the screening masses), it has 
been exploited in the method of "dimensional reduction" (see [41], 4^, 43, 0] and references 
therein). We emphasize, however, that the damping rate is a dynamical quantity, and 
the usefulness of the dimensional reduction for this problem is not a priori obvious, given 
the subtleties of the analytic continuation from Matsubara to real external energy. If 
a dimensional reduction occurs in the computation of the large time behaviour, this is 
because of the particular IR behaviour of the magnetic photon propagator, as displayed in 
eqs. (|2.16|) or ( 2.29|) . The dynamical information which is contained in the later equations 
refers not only to the absence of the magnetic screening, but also to the phenomenon of 
Landau damping. 

It is also worth emphasizing that our result takes into account only the most singular 
terms of the perturbative expansion. Because of the approximation used, we have lost 
control on the subleading terms. Although, in a strict perturbative sense, these are a 
priori less important, one cannot completely exclude the possibility that they may still 
modify our results in a qualitative way. It is hard to see however how they could destroy 
the quasiparticle picture, which we have shown to survive after a complete treatement 
of the leading IR divergences. Improvements of our solution may require an appropriate 
generalisation of the Bloch-Nordsieck model to finite temperature, a task that we have 
explored in this paper, but without reaching a definite conclusion. There are at least two 
points where the thermal BN model could possibly complete our previous analysis: the 
dynamical emergence of the upper cut-off ~ gT (recall that, in the effective 3-dimensional 
theory, this cut-off has been introduced by hand), and, related to this, the consistent 
computation of the subleading terms beyond \n{upt) in eq. ( ^.13| ), that is, the terms of 
order g'^T which multiply the time in the exponent of A(t). 

It is finally natural to ask what is the relevance of the present solution for the 
non-Abelian QCD plasma. It is widely believed that the self-interactions of the chro- 
momagnetic gluons may generate magnetic screening at the scale g'^T. As a crude 
model, we may include a screening mass /x ~ g^T in the magnetostatic propagator in 
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the QED calculation. Then eq. ( |3.42| ) provides, at very large times t ^ 1/g'^T, an ex- 
ponential decay, A^(t) ~ exp(— 7^) with 7 = aTln{up/ fi) = aT\ia{l/g). However, in 
the physically more interesting regime of intermediate times l/gT ^ t ^ 1/g'^T, the 
behavior is governed uniquely by the plasma frequency, according to our result ( |4.13| ): 
A^(t) ~ exp{—aTt In ujpt). Thus, at least within this limited model, which is QED with 
a "magnetic mass", the time behavior in the physical regime remains controlled by the 
Bloch-Nordsieck mechanism. But, of course, this result gives no serious indication about 
the real situation in QCD, since it is unknown whether, in the present problem, the effects 
of the gluon self-interactions can be simply summarized in terms of a magnetic mass. 
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Appendix A 

In this Appendix, we collect the sum rules for the photon spectral densities which 
are used in section 2.2. 

The electric and magnetic spectral densities are defined in eq. ( p.l9| ) in terms of 
the corresponding propagators. In the hard thermal loop approximation, they involve 
both pole and cut pieces, as shown in eq. ( |2.22| ). They satisfy the following sum-rules P^, 
which trade the integrals over the off-shell spectral densities Pi,t{(lo,(l) for functions of 
ujs{q) and Zs{q) : 

dqo 1 1 zi{q) 

Awo,g) - 



g Svrgo ' q'^ + m% ujf{q)' 

:r—Pt[qo,q) - - 
-q ZTcqo q Uf [q) 

^qoPt{qo,q) = 1 - zt{q). (A.l) 
-q zn 

The first two of these sum rules are obtained by simply setting = in the spectral rep- 
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resentations ( p.l8| ), and by using *Ai{0,q) = — + m|)), *At(0,g) = — together 
with eq. (|2.22| ). As for the third one, this is obtained by inserting eq. ( |2.22| ) into the 
famihar sum-rule 



/ -^Qo pt{qo,q> = 1, 



(A.2) 



which is a consequence of the equal-time commutation relation for the quantum fields 13 



The use of the sum rules ( [A.l|) is convenient to study both the ultraviolet and 
the infrared behavior of the g-integral in eq. (|2.26|) . To this aim, we need the dispersion 
relations uji^tiq) IH) 0; and the corresponding residues zi^t{q), which, in our conventions, 
read: 



2^2 (^2 _ ^2 



3a;2^2 _ (^2 _ ^2)2 



2u;f{ujf/q' - 1) 

3.;2 _ (^2 _ g2) 



At large momenta, q 3> ujp, we have the approximate expressions [SC 



(A.3) 



Zt{q) ^ 1 



H 

4g2 



In 



3^2 



-3 



uj^i{q)-q\l + Axi{q)), 
8g2 



where 



xi[q) = exp I -T-^ - 2 



3^2 



From eqs. ( |A.1| ) and (|A.4|) , we obtain, for q':^ ojp (recall that m|, = 3a;p), 



-q 2nqo 



-q 27rgo V Q J 



2g4 



(A.4) 



(A.5) 



(A.6) 



These estimates show that the integrand in eq. (|2.26|) behaves like oj'^/q^ for momenta 
q > ujp. 

We turn now to momenta q <^ ujp. We then have 



u;2(g) ~u;^ + 6g75, 

-,2 



1 - 



9 

5^2 



(A.7) 
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so that 



1 dgo 



1 dgo 



A(go,g) 
A(go,g) 



/ v^goA(go,g) 



4 1 

15 ^' 



(A.8) 



When these expressions are inserted in eq. ( |2.26| ), the contribution in of the magnetic 
spectral function (the second hne in eq. (|A.8|)) generates a logarithmic IR singularity. 



Appendix B 

Since there is no phase-space available for the direct decay of the on-shell fermion 
into a pair of massless particles, one expects that the damping rate computed from the 
bare one-loop fermion self-energy should vanish. However, at finite temperature, this 
argument is complicated by infrared singularities which arise because of the enhancement 
of collinear singularities by the Bose-Einstein thermal factor. 

To illustrate this problem, we consider the calculation of the damping rate to 
bare one-loop order in the Coulomb gauge. This is obtained by simply replacing, in 
eq. ( 2.21|) , the photon spectral functions with their bare counterparts, namely p^^^ = 



and pf'\qo,q) = Po{qQ^q), with po from eq. ( |1 . 2| ) . In the on-shell limit, the whole contri- 
bution to 7 comes from space-like photons, with |go| < q- However, since the free spectral 
density (|1.2| ) has support precisely at the integration limits go = we should be more 
careful when evaluating eq. ( f^.21D in the on-shell limit uj ^ p. For uo close to, but different 
from, p, the latter equation yields (compare with eq. ( p.25|) ) 



p) ^tt/T f ^ r 6{uj-p-q, + qcose){l-cos'e)po{qo,q). {B.l] 



After the angular integration, we obtain 

g'^T /-"^ dg r^-p+q dqo 



g I f°° dq f^-p+q aqo r 
7o(^~p) ^ ——{uj-p) — {2qo - {uj - p)) 6{qo - q) - 6{qo + q) 

OTT JO q^ Ju-p-q qo 
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If we let now u p, the factor in the front of the last integral goes to zero, but the integral 
itself becomes IR divergent. An explicit calculation shows that the r.h.s. of eq. ( |B.2| ) is 
in fact independent of {uj — p), and equals 

70 = ^. (B.3) 

This result is, however, unphysical. It arises from the emission or the absorbtion of 
collinear {6 = 0, Qq = q or, respectively, 6 = n, qo = —q) massless photons, whose 
contributions are enhanced by the Bose-Einstein factor T/go- Such contributions do not 
survive screening corrections. However, since the gauge-dependent terms in the photon 
propagator are not modified by the plasma effects, they may generate — by the mechanism 
alluded to before — a non- vanishing contribution to the on-shell self-energy Note 
that an entirely similar problem arises in the three-dimensional gauge theories at zero- 



temperature, when computing the dispersion equation to one loop order] 45 



To overcome this problem, it has been suggested p5| , ^ to take the on-shell limit 
in the presence of an IR regulator, say, an IR cut-off /i. With such a cut-off, the g-integral 
in eq. ( p.2|) remains finite as uj p, and the total result for joluj = p) vanishes. Thus, 
the damping rate remains zero at the bare one-loop level, as expected. In the same way 
one verifies that the dispersion relation is gauge- independent, as it should [p9|]. On the 
other hand, the residue of the propagator at the mass-shell becomes dependent on the IR 
cut-off /i, and (linearly) divergent as — 0. 



Appendix C 

We verify here, on an explicit two-loop calculation, some general features of the 
infrared behavior of the on-shell self-energy in perturbation theory. Specifically, we shall 
show that the leading divergences are power-like, and can be fully taken into account by 
restricting all the internal Matsubara sums to their zero frequency photon modes. 

At two loop order, the fermion self-energy is given by the two diagrams in Fig. 0, 
which yield 

[ixSoip + qhu + luSoip + k)^x] "D^^iq) *D''\k) . (C.l) 
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Figure 11: Two-loop diagrams for the fermion self-energy. 



The notations here are as in eq. (|2.17| ); for instance, p° = iUn = i{2n + l)7rT, = ium = 
i27rmT and = iur = i27TrT, with integers n, m and r. According to eq. ( p.. 61 ), the correc- 
tion to the positive mass-shell is determined by the function Tj^{uj,p) = tr (ji-i-{'p)Tj{uj , p))/2. 
As C(j ~ the most singular contributions to are obtained by using the effective Feyn- 
man rules described at the end of section 3.1. At two loop level, this amounts to replacing 
eq. dOl]) by 



(27r)3y (27r 
[Go{p + q) + Go{p + k)]D{q)D{k) 



Goip + q) Goip + q + k) 



(C.2) 



where D{q) = v'' * D^^ {q)v^ is the (HTL resummed) propagator of the magnetic photon. 
(We recall that the electric propagator does not yield IR singularities.) Eq. (|C.2|) is 
precisely the two-loop self-energy in the Bloch-Nordsieck approximation. 

The Matsubara sums over Um and uJr are conveniently performed by contour meth- 
ods, and by using the spectral representation ( p.l8| ) of *At{q). In doing this, one gets sev- 
eral terms, corresponding to the poles of the various propagators in the complex planes 
and k^. Every such a term involves three energy denominators, and the product of two sta- 
tistical factors. The latter are either of the bosonic or of the fermionic type, according to 
whether they correspond to poles of a photon propagator, or of an electron propagator, re- 
spectively. When the external energy approaches the tree-level mass-shell, uj p = \ ■ p, 
all the energy denominators are soft, of the type l/(go — v ■ q), and may give infrared 
problems. (The hard energy denominators, which were potentially present in the full two- 
loop self-energy (|C.1|) , have been eliminated by the simplified Feynman rules leading to 
eq. (|C.2| ).) Then, the leading IR singularities arise uniquely from the terms which involve 
the product of two Bose-Einstein distribution functions, since N{q^)N{k^) ~ T"^ /{q^k^) at 
soft momenta. By isolating these most singular terms, we obtain, after a straightforward 



51 



calculation, 



d^q f d^k dqo 
(2^ J (27r)3 i_oo 2^ 
1 



27rA;n 



*Ptiko,k) 



1 



+ go - V ■ (p + q) 
1 



go + fco - V ■ (p + q + k) 
1 



_|_ I tQ 3^ 

_co' + go - V ■ (p + q) + A;o - V ■ (p + k) _ 
where it is understood that the external energy carries a small positive imaginary part 
uj + irj) . 

The energy integrals over go and over /cq involve both the pole and the cut pieces of 
the photon spectral density. However, it is only the off-shell (or cut) piece of *pt which 
yields a singular contribution, so we may as well restrict the aforementioned energy inte- 
grals to space-like momenta, |go| < g and |A;o| ^ and replace the full spectral functions 
by Pf Then, the subsequent analysis follows closely the discusion of the (resummed) one- 
loop self-energy in section 2.2. The singular domain is that of very soft photon momenta, 
q, k gT, where we can use eq. ( p^.29| ) to replace Ptilo ^ <l)/<lo by (27r/g^)5(go). At the 
same time, we have to supplement the momentum integrations with an upper cut-off of 
gT. The net effect is that the leading singular piece of ^^^\uj ~ p) is 



the order of oOp ~ 



the same as it would be obtained by retaining only the static terms Um = 
Matsubara sums of eq. (|C.2| ). That is. 



in the 



^(2) 



p) 



d^q 1 f d^k 1 



(27r)3 g2 
1 



1 



V ■ (P + qj 



(27r)3 P u 
1 

uj — \ ■ 



ip + q) 



ip + q 



+ 



k) 

1 



(P + k) 



(C.4) 



Since this is divergent as u 
IR cut-off fi, and obtain 



V ■ p, we take the mass-shell limit in the presence on an 



d^q 1 



(27r)3 g2 (v ■ q — zr/)^ J (27r)3 A;2 v ■ k — ir] 



d^A; 1 



. 2 {aTf ujp 
I In — 



(C 



We thus find the linear plus logarithmic infrared divergence mentioned in section 2.3. 

According to eq. ( p.33| ), the computation of 7*^^^ — the two-loop contribution to the 
damping rate — requires also the one-loop residue, z^'^^{p) — 1 = {dTj^_^^ / dun) . Similarly 
to eq. ( p.4|) , we obtain the leading IR-singular contribution to S^^"* in the form 



d^q 1 



(27r)3 g2 00 
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V ■ (p + q) 



17] 



(C.6) 



and thus 

^(1) p _ 1 ^ _g^T / — ^ - — ^ , C.7 

in the presence of the IR regulator. The hnear IR divergence of the residue compensates 
the dominant singularity of the two-loop self-energy ( |C.5| ), so that the leading contribution 
to 7*^^-* — which remains beyond the accuracy of the present computation — is of the order 
{a^T"^ /ujp){\n.{ijjp/ ~ 9^T{\n.{ijjp/ ^))'^ . Even if still divergent as — > 0, this does not 
contribute to the order g'^T which is our concern here. 

Let us finally provide an all order argument for the cancellation of the strongest, 
power-like, infrared divergences in the perturbative evaluation of 7. To this aim, we 
consider the Dyson-Schwinger equation for the fermion self-energy within the effective 
three-dimensional Bloch-Nordsieck theory: 

p) = -/T / ^ 5(0;, p + q) F(p + q, p) Z}^^(q), (C.8) 
where D]^{(\) = S is the full BN propagator, 

5(^, p + q) = f ' 

- V ■ (p + q) - S(cj, p + q) 
and P (p + q, p) is the full vertex, which is related to S via the Ward identity 

q^T^{p + q,p) = S-\uj,p + q)-S-\uj,p). (C. 10) 



We make now the usual assumption ||53[ that the dominant IR behavior involves only the 



longitudinal piece of the vertex. This is entirely determined by the Ward identity: 

r^'(p + q,p) = ^\S-\uj,p + q)-S-\iU,p)]. (C.ll) 
V ■ q L J 

When inserted in eq. (|C.8|) , this yields 

d^q 1 

{27c)^ q^ V ■ q — it] 

As already explained, eq. ( P.8|) reproduces the most singular terms of the perturbative 
expansion, and this remains true after inserting the approximation ( |C.11| ) for the vertex 
function, as can be verified explicitly by developing eq. ( |C.12| ) in perturbation theory. We 



J:{u,p) ^ -g'T /t^^ ^—[^ - S-\u,p)S{u,p + cD]. (C.12) 



now take the on-shell limit in the presence of an IR cut-off /i, taken as a small photon mass. 
As long as /i 7^ 0, there is no IR problem, and we expect the mass-shell to correspond 
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to a simple pole of the exact propagator. Thus, S~^{u!,p) vanishes on shell, and the 
second term in eq. (|C.12D gives no contribution. The leading contribution to the on-shell 
self-energy reads then 

E(on - shell) ^ -^T / -fi^ \ . ^ ^ ^ In , (C.13) 

and coincides with the IR singular part of the one-loop self-energy. This is only possible if 
the aforementioned compensation of the leading power-like divergences holds in all orders. 
Note that the above arguments become meaningless in the physical limit fi 0, where 
not only does the estimate ( |C.13|) become logarithmically divergent, but the integral 
multiplying S'~^(ci;,p) also diverges on the mass-shell. 
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